Local scale-invariance and ageing in noisy systems 



Alan Picone and Malte Henkel 

Laboratoire de Physique des Materiaux/ Universite Henri Poincare Nancy I, 
B.P. 239, F - 54506 Vandceuvre les Nancy Cedex, France 

Abstract 

The influence of the noise on the long-time ageing dynamics of a quenched ferromagnetic spin 
system with a non-conserved order parameter and described through a Langevin equation with a 
thermal noise term and a disordered initial state is studied. If the noiseless part of the system is 
Galilei-invariant and scale-invariant with dynamical exponent z = 2, the two-time linear response 
function is independent of the noise and therefore has exactly the form predicted from the local 
scale-invariance of the noiseless part. The two-time correlation function is exactly given in terms of 
certain noiseless three- and four-point response functions. An explicit scaling form of the two-time 
autocorrelation function follows. For disordered initial states, local scale-invariance is sufficient for 
the equality of the autocorrelation and autoresponse exponents in phase-ordering kinetics. The 
results for the scaling functions are confirmed through tests in the kinetic spherical model, the 
spin-wave approximation of the XY model, the critical voter model and the free random walk. 
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1 Introduction 



The study of ageing phenomena as they are known to occur in glassy and non-glassy systems presents 
one of the great challenges in current research into strongly coupled many-body systems far from 
thermal equilibrium. A common example of this kind of system is obtained as follows. Consider a 
magnet at a high-temperature initial state before quenching it to a final temperature T at or below 
its critical temperature Tc > 0. Then the temporal evolution of the system with T fixed is studied. 
A key insight has been the observation that many of the apparently erratic and history-dependent 
properties of such systems can be organized in terms of a simple scaling picture jHOl- Underlying 
this phenomenological picture is the idea that the ageing phenomenon and the related slow evolution 
of the macroscopic observables comes from the slow motion of the domain walls which separate the 
competing correlated clusters. The domains are of a typical time-dependent size with length-scale L{t), 
see El m inn I2ni for reviews. In recent years, much work has been performed on the ageing 
phenomena of simple ferromagnetic systems, in the hope that these systems might offer insight useful 
for the refined study also of ageing glassy materials. It has turned out that ageing is more fully revealed 
in two-time observables, such as the two-time (auto-)correlation function C{t, s) or the two-time linear 
(auto-)response function R{t, s) defined as 

5{m 



C{t,s):=m)<P{s)) , R{t, 



6h{s) 



(1.1) 

h=0 



where denotes the time-dependent order parameter, h{s) is the time-dependent conjugate magnetic 
field, t is referred to as observation time and s as waiting time. One says that the system undergoes 
ageing if C or i? depend on both t and s and not merely on the difference t = t — s. According to the 
dynamical scaling alluded to above, one expects for times t, s 3> tmicro and t — s ^ ^micro, where tmicro 
is some microscopic time scale, the following scaling forms 

Cit, s) = s-'fcit/s) , Rit, s) = s-'-'^fnit/s) (1.2) 

such that the scaling functions fc,R{y) satisfy the following asymptotic behaviour 

fc{y) ~ y-'^/^ , My) ~ y''''^' (i-S) 

as y — ^ oo and where Ac and Xr, respectively, are known as the autocorrelation [2711^11 and autoresponse 
exponents j77j, a and b are further non-equilibrium exponents and z is the dynamical exponent, where it 
has been tacitly assumed that the typical cluster size grows for late times as L(t) ~ t^/^. The derivation 
of such growth laws from dynamical scaling has been studied in great detail [85J. For a non-conserved 
order parameter and T < Tc, the dynamical exponent z = 2. The exponents Xc,r are independent 
of the equilibrium exponents and of z jSHl IHH IIH]- Since a long time, the equality Ac = Xr had 
been taken for granted but recently, examples to the contrary have been found for either long-ranged 
initial correlations in ageing ferromagnets |77j or else in the random-phase sine-Gordon (Cardy-Ostlund) 
model [89j. On the other hand, a second-order perturbative analysis of the time-dependent non-linear 
Ginzburg-Landau equation reproduces Ac = A^ ^7\. The precise relationship between Ac and Xr 
remains to be understood. If one uses an infinite-temperature initial state, one has Xc = Xr> d/2 [93] . 

For ageing ferromagnetic systems with a non-conserved order-parameter, the value of the exponent 
a depends on the properties of the equilibrium system as follows jl71|S2|. A system is said to be in class 
S if its order-parameter correlator Ccq{r) ~ exp(— |r|/^) with a finite ^ and it is said to be in class L if 
Ceq(r) ~ |r|-('^-2+^), where ?7 is a standard equilibrium critical exponent. Then 

1/z ; for class S ,^ ^, 

{d — 2 + ri)/z ; for class L 
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For example, in c? > 1 dimensions, the kinetic Ising model with Glauber dynamics is in class S for 
temperatures T < T^. and in class L at the critical temperature T = T^. It is generally accepted that 
6 = for T < Tc and b = a if T = Tc, see e.g. [39j. 

The distance from equilibrium is conveniently measured through the fluctuation- dissipation ratio 

iiniini 

X(t,5):=Ti?(t,.)(^^^^) ' (1.5) 

At equilibrium, the fluctuation-dissipation theorem states that X{t, s) = 1. Ageing systems may also 
be characterized through the limit fluctuation-dissipation ratio 

Xoo = lim (\im X{t,s)) (1.6) 

Below criticality, one expects X^o = 0, but at T = the value of X^o should be universal according 
to the Godreche-Luck conjecture jHTl I38j. This universality has been confirmed in a large variety of 
systems in one and two space dimensions [SHI El 1^ EH]- The order of the limits is important, since 
limt^oo {lim.s^ooX{t, s)) = 1 always. 

While these statements exhaust the content of dynamical scaling, it may be asked whether the 
form of the scaling functions fc,R{y) might be fixed in a generic, model- independent way through a 
generalization of that symmetry. Indeed, it has been shown that an infinitesimal global scale- 
transformation t \—>- {1 + e)^t, r {1 -\- e)r with a constant e such that \e\ <^ 1 can for any given 
value of z be extended to an infinitesimal local scale-transformation where now e = 6{t, r) may depend 
on both time and space. ^ It can be shown that the local scale-transformations so constructed act 
as dynamical symmetries of certain linear field equations which might be viewed as some effective 
renormalized equation of motion. Practically more important, assuming that the response functions 
of the theory transform covariantly under local scale-transformations, the exact form of the scaling 
function fR^y) is found jlHlllS] 

UV) = rov'^^-'^^Hv ~ ly'-" (1.7) 

where tq is a normalization constant.^ Indeed, in deriving this result one actually only requests that 
R transforms covariantly under under the sttfealgebra of the infinitesimal local scale-transformation 
which excluded time-translations. We say that a theory where the n-point functions built from certain 
'quasiprimary' fields transform covariantly under an algebra of such extended scale-transformations 
is locally scale-invariant jJUl |33] . The prediction (jl.7|) has been confirmed in a large class of ageing 
ferromagnets as reviewed in j^UlEOl. The status of the scaling function fc{y) of the spin-spin correlator 
is less clear, however. Building on the Ohta-Jasnow-Kawasaki approximation (see [H]) Gaussian closure 
procedures jHllH3| in the 0(n)-model produce approximate forms for fc{y) but we do not know of any 
other approach which does not involve some uncontrolled approximation. 

Given the phenomenological success of ()1.7p . we wish to understand better where such a supposedly 
general and exact result derived from a dynamical symmetry and without using any model-specific 
properties might come from. In this paper, we shall concentrate on the important special case z = 2 
which describes the phase-ordering kinetics after a quench to a temperature T < Tc of a ferromagnet 
with a non-conserved order-parameter |85]. We recall that local scale-transformations are dynamical 

^This extension of dynamical scaling has an analogue in critical equilibrium systems: there global scale invariance can 
be extended to conformal invariance, see ,17, .24, 44 for introductions. 

■^In order to avoid misunderstandings, we recall that 1)1. 7|l holds for the total response function as defined in 1)1. 1|) 
without any subtractions meant to extract an 'ageing part'. 
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symmetries of certain differential equations, such as the free diffusion/free Schrodinger equation for 
z = 2. Indeed, the maximal dynamical symmetry of the free diffusion equation dtcf) = DA(f) is known 
since a long time to be the so-called Schrodinger group jHSl HH EO], to be defined in section 2. Also, 
it is well-established that the same group also describes the dynamic symmetry of non-relativistic free- 
field theory |121EZI- It also arises as dynamical symmetry in certain non-linear Schrodinger equations 
ISH IS21 ins ESI EI] , the Burgers equation [^EEI or the equations of fluid dynamics [Hj. If in addition D 
is also considered as a variable, Schrodinger invariance in d dimension becomes a conformal invariance 
in (i+2 dimensions jlH^^. The classification of non-linear equations and of systems of equations admitting 
as a dynamical symmetry the Schrodinger group or one of its subgroups (e.g. the Galilei group) has 
received a lot of mathematical attention, see jHOl E21 EHl EH E21 UHl IZ2j- The extension to dynamical 
exponents z 2 needed for quenches to criticality or for glassy systems will be left for future work. 

However, the setting just outlined is not yet sufficient for the description of ageing phenomena. 
Rather, we are interested in the time-dependent behaviour of spin systems coupled to a heat bath at 
temperature T. It is usually admitted that after coarse-graining, this may be modeled in terms of 
a Langevin equation. If there are no macroscopic conservation laws, the Langevin equation for the 
coarse-grained order parameter = (f){t, r) should be model A in the Hohenberg-Halperin classification 

m 

where 7i is the classical Hamiltonian, and D stands for the diffusion constant or equivalently some relax- 
ation rate. Thermal noise is described by a Gaussian random force t] = r]{t, r) and is thus characterized 
by its first two moments 

{ri(t,r)) = 0, {ri{t,r)ri{s,r')) = 2DT5{t~ s)5{r -r') (1.9) 

where T is the bath temperature. It is well-known [SSI El that this formalism describes the relaxation of 
the system towards its equilibrium state given by the probability distribution Peq ~ e~^/^. In addition, 
initial conditions must be taken into account and are described in terms of the initial correlation function 

a(r - r') := C(0, 0; r, r') := (0(0, r)0(O, r')) (1.10) 

and where we already anticipated spatial translation invariance. 

Neither the thermal noise nor the initial correlations described by a{r) are included into the local 
scale-transformations as studied in |46] which come from systems such as the free diffusion equation. In 
this paper, we want to show how both these sources of fluctuations may be taken into account and we 
shall explicitly derive the two-time response and correlation functions. Our analysis will be restricted 
to the case where z = 2 which for instance is already enough to describe ageing below criticality. 

As we shall show, it is useful to slightly generalize the problem and to consider the kinetics of systems 
which in the simplest case may be described by a quadratic Hamiltonian of the form 



n[<P] = ^J dtdr 



dr 



;i.ii) 



where v{t) is a time-dependent external potential. Formally, at the level of relativistic free- field theory, 
v{t) corresponds to a (time-dependent) mass squared which would measure the distance from a critical 
point. Alternatively v{t) may be viewed as a Lagrange multiplier in order to ensure the constraint 
(0(t, r)(j){t, r)) = 1 and we shall make this explicit through the example for the kinetic spherical model 
in section 5. In a physically more appealing way, time- dependent potentials arise when a many-body 
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system is brought into contact with a heat bath whose temperature T{t) is time-dependent In 
this paper, we shall be interested in the dynamics symmetries of Langevin equations derived from a 
free- field Hamiltonian In particular, we shall compare the situation without (i.e. T = 0) and 

with thermal noise (i.e. T > 0). For simplicity, we shall refer to all equations ()1.8|1 obtained from the 
Hamiltonian (jl.llj) as 'free Schrodinger equations'. 

The study of the dynamic symmetries of such free-field theories will yield useful insights which we 
expect to extend to physically more realistic interacting field theories where Ti would also contain higher 
than merely quadratic terms. If we identify = 2im, it is clear that the order parameter (j) is given 
by a noisy Schrodinger equation in an external time-dependent potential v{t). 

This paper is organized as follows. In section 2, we first review the basics of Schrodinger-invariance 
in the absence of thermal noise and without initial correlations and then show that through a gauge- 
transformation involving v{t), the entire phenomenology of ageing and in particular ()1.7|) can be re- 
produced. We also consider the selection rules which follow from Galilei-invariance. In section 3, 
after having reformulated the problem in terms of the field-theoretic Martin-Siggia-Rose formalism, we 
study the effects of thermal noise and/or initial correlations on free-field theory given by a Hamiltonian 
In section 4, these results are extended to any field theory with for T = and a(r) = is 
Galilei-invariant. We find that the two-time response function R is independent of both T and a(r) 
and obtain a new reduction formula ()4.9|) which relates C to certain three- and four-point response 
functions to be evaluated in the noiseless theory and discuss the scaling of the resulting two-time auto- 
correlation function. In sections 5-7, these results are tested in several exactly solvable systems (with 
an underlying free-field theory) undergoing ageing with z = 2, namely the kinetic spherical model, the 
XY-model in spin-wave approximation, the critical voter model and the free random walk. We conclude 
in section 8. Appendix A deals with technical aspects of Gaussian integration and appendix B analyses 
a special four-point response function. In appendix C we consider a generalized realization of local 
scale-invariance and its application to the ID Glauber-Ising model 



2 Local scale-invariance: a reminder 
2.1 Schrodinger-invariance 

We begin by reviewing the kinematic symmetries of Schrodinger equations with a time-dependent 
potential, but without a noise term. A long time ago, Niederer obtained the maximal kinematic 
symmetry group of the Schrodinger equation for an arbitrary potential v = v{t,r) and he also gave 
a few examples where that group is isomorphic to the maximal kinematic group Sch[d) of the free 
Schrodinger equation inHUZOl- The group Sch{d) is called the Schrodinger group (TUj . 

We recall the definition of Sch{d). On the time and space coordinates (t, r) it acts as {t,r) 
{f, r') = g{t, r) where 

, at + (3 , Ur + vt + a 
t\ — >t = ——— , ri — >r = — — ; a6 - = 1 (2.1) 

•yt + d 'jt + 

where 7?. is a rotation matrix. The action of Sch{d) on the space of solutions (p of the free Schrodinger 
equation is projective, that is, the wave function (p = (f){t, r) transforms into 

0(t, r) (T,0) (t, r) = f,[g-\t, rMg~\t, r)] (2.2) 

and the companion function fg is explicitly known jJOl EI]- The projective unitary irreducible rep- 
resentations of Sch{d) are classified [ZOj- We now carry this over to field theory and consider fields 
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transforming according to (|2.2|) . By analogy with an analogous terminology in conformal field theory 
[3j, a field (f) transforming according to ()2.2j) and with g{t,r) given by ()2.H) is called quasiprimary |46j. 
Schrodinger-invariance is the z = 2 special case of local scale-invariance, with time-translations added. 

Besides the examples given in [7T] , there exist further noiseless Schrodinger equations with a maximal 
kinematic group isomorphic to Sch{d). Consider the noiseless Langevin equation 

D^i^^^llTl = A0(t, r) - vimt, r) (2.3) 

where D is the diffusion constant. This equation can be reduced to the free Schrodinger equation 
D^^dt^{t,r) = A\E'(t, r) through the gauge transformation 

(j)[t, r) = ^(t, r) exp (^-D duv{u)^ (2.4) 

Since the kinematic symmetries of the free Schrodinger equation are well-understood and the realization 
of the Schrodinger group is explicitly known, the corresponding realization for the case at hand, similar 
to ()2.2|) . readily follows. 

It turns out that the only change occurs in the companion function fg. Let ff'^ stand for the 
companion function of the free Schrodinger equation, then because of the gauge transformation eq. ()2.4j) 
we find 

/,(t,r) = /f (t',r)exp ^ duv{u)^ (2.5) 

where t' = t'{t) has been defined in eq. ()2.H) . The generators of the Lie algebra scf)^ of this realization 
of Scli(l), appropriate for the equation (j2.3|) with the potential f (t), read 

X_^ = -dt + ^ time drift 
2M 

1 X vit] 

Xn = —tdt rdr h t — — dilatation 

" 2 2 2A^ 

Xi = —t^dt — trdr — xt + ^^TTT-T — -^^^ special Schrodinger transformation 

Yli/2 = —dr space translation 

Yi/2 = —tdr — Air Galilei transformation 

Mo = -M phase shift. (2.6) 

where we expressed the diffusion constant D^^ = 2A4 as the 'mass' Ai and x denotes the scaling 
dimension of the wave function (pit^r). Of course, x = d/2 for a solution of the free Schrodinger 
equation, but it will be useful to consider arbitrary values of x as well. 

The non- vanishing commutators of the Lie algebra sc^i spanned by the generators ()2.(ij) are 

[X„, X„,] = {n- n')X„+„, , [X„, F„] = - m) , [F1/2, ^-1/2] = Mq (2.7) 

where n, n' G {±1, 0}, m G {±|} and with straightforward extensions to d > 1, see jl^ . 

2.2 Galilei-covariance of correlators 

When discussing the dynamic symmetries of a time- dependent statistical system, the requirement of 
Galilei-invariance plays a particular role. Indeed, for a system with local interactions and which is 
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invariant under space translations, scale transformations with a dynamical exponent z = 2 and in 
addition Galilei-invariant, it can be shown that there exists a Ward identity such that the system is 
also invariant under the 'special' Schrodinger transformation ()2.6p [121113 EH]- 

We shall be particularly interested in the two-point correlator C and the linear response function R 
built from the order parameter (j){t,r). Using Martin-Siggia-Rose theory (MSR theory) which we shall 
briefly review in section 3, these may be expressed in terms of and the so-called response field as 
follows 

C{t,s;r,r') := m,r)4>{s,r')) 



= (0(t,r)0(s,r')) (2.^ 

h=0 



R{t,s;r,r') : 

oh{s, r' 

C{t,s-r,r') := {'^{t,r)^{s,r')) 

where h is the magnetic field conjugate to the order parameter 0. Later, we shall often refer to t as the 
observation time and to s as the waiting time. 

Generalizing the above definition to v{t) ^ 0, we say that a field is quasiprimary if its infinitesimal 
change under scf)^ is given by the generators ()2.6p with v(t) 0. A quasiprimary field is characterized 
by its scaling dimension x and its 'mass' A4 > 0. In turn, if the response field associated to is also 
quasiprimary, it has a scaling dimension denoted by x and the 'mass' 

M = -M<0 (2.9) 

This important fact will be used later on. The argument leading to the result (j2.9p was discussed in 
detail in [48^ and will not be repeated here. 

If both and transform as quasiprimary fields of a Schrodinger-invariant theory, the generators 
eqs. ()2.6|) can be used to derive restrictions on the form of any multipoint correlator and in particular 
determine the two-point functions completely. If Xi is any of the generators of 5cf)^ acting on the i^^ 
particle in a n-point correlator 2t{tj,rj} where i = 1, . . . ,n (see (j2.8|) for n = 2), we have a set of 
differential equations 

+ ... + Xn)'^{ti,ri} = (2.10) 

If rotation invariance can be assumed (and we shall implicitly do so throughout this paper), for the 
calculation of the two-point functions it is enough to consider the one- dimensional case and use the 
generators of eq. (|2.6|) . Then a straightforward calculation jSU] gives, provided the 'mass' of the order 
parameter is positive A4 > 0, see e.g. [ISl liU] 

Coit,s;r,r') = 

Coit,s;r,r') = (2.11) 

whereas the response function is basically the gauge-transformed expression of the well-known zero- 
potential Gaussian response TZ, i.e. 



R^{t,s;r,r') = '^n{t,s;r,r') 
k{s) 

TZ (t, s; r, r') = 5x,x roQ{t~ s) {t - 3^ exp 

1 



Mir- r'? 



2 {t 

kit) := exp(-^ / duv{u)] (2.12) 
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where tq is a normalization constant and G is the Heaviside function which expresses causahty. As they 
stand, eqs. ()2.11I2.12|) hold for T = and we shall from now on use the index to remind the reader 
of this fact. 

On the other hand, if the system is not rotation-invariant, we can repeat the same argument in 
any fixed direction of space and the non-universal constant Ai becomes direction-dependent. Indeed, 
rotation invariance is broken for phase-ordering systems defined on a lattice [SniEZl for sufficiently small 
temperatures. Even then, local scale invariance still holds in every single space direction, as exemplified 
in the 2D and 3D Ising models with Glauber dynamics j49j . 

A few comments are in order. 

1. Eqs. ()2.11I2.12|) provide a manifest example of the superselection rule of Galilei invariance, also 
known as Bargman superselection rule |2]. Explicitly, if r^) are Galilei-covariant fields, each 
with a 'mass' Aii, Galilei-covariance implies [21 ESI 

(<l>i(ti, ri) . . . <l>„(t„, r„)) = 6M^+...+M„fi F {{U, r^}) (2.13) 

By physical convention, the 'masses' of the fields are non-negative, viz. }Ai > 0. Furthermore, 
the response fields have negative 'masses' Aii < and the result 1)2.111) follows. 

2. In the introduction, we reviewed the result p.7p for the autoresponse function, derived for arbi- 
trary z from a generalization of Schrodinger invariance jlHl ESI- For z = 2, eq. ()1.7|) coincides 
with our result ()2.12j) provided that 

X = X = 1 + a 
vit) = i2M) ^ + (2.14) 

3. However, there is an important difference in our derivation of the scaling form of Ro{t, s) with 
respect to [lUJ : because time-translation invariance is broken in ageing phenomena, covariance 
of Rq was required to hold merely under a subalgebra of scf)^ where the time-translations were 
excluded. Indeed, in jlH], such subalgebras were studied systematically and a relationship with 
the parabolic subalgebras of the (complexified) conformal algebra conf^_^_2 was found. On the other 
hand, the realization of sci)^ used in j45ll3SlllHl applies to a vanishing potential v(t) = 0. 

Here, we do not follow that point of view. We consider the more general realization of the entire 
algebra sci)i with a time-dependent potential v{t) and require that both (p and (p transform as 
quasiprimary fields under the whole set of generators ()2.6|) . 

That these two different approaches, given the conditions (|2.14j) . yield the same phenomenology 
of the scaling of the two-time autoresponse function is our first result and will be crucial for the 
developments to follow. 

4. At first sight, the result (|2.1ip that the two-time autocorrelator C(t, s) = {4>(t, t)0(s, r)) = may 
appear strange and indeed does not hold true in concrete models. In the next sections, we shall 
show that this apparent contradiction comes from the fact that the noiseless Schrodinger equation 
does not take the thermal noise into account. As we shall see, the reformulation of Schrodinger 
invariance in ageing systems in terms of a noisy Schrodinger equation with a time-dependent 
potential allows to arrive at physically meaningful predictions for correlation functions. Explicit 
confirmations in exactly soluble models will be presented. 
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3 Response and correlation functions for non-interacting Gaus- 
sian theories 



3.1 The Martin-Siggia-Rose formalism 

It is useful to treat noisy Langevin equations in the context of the Martin-Siggia-Rose (MSR) formahsm 
jHHl Eni inn 1211 m] • in equihbrium, the integration over the Gaussian noise rj can be carried out by 
introducing a response field 0.^ It can be shown that the stochastic Langevin equation ()1.8|1 can be 
obtained from the following effective action S[0,(/)] 

E[0,0] = jAtAr4>{^ + D^-^^-\ ldtdrdt'dr'{v{t,rUt',r'))^it,r)^it',r') (3.1) 

This action appears in the generating functional Z = J VcjyDcj) e"^!*^'*^' expressed as a path integral. 
In this way, the original dynamical problem in d dimensions has been mapped onto one of statistical 
mechanics m. d+1 dimensions. 

As long as one is merely interested in equilibrium behaviour and provided the dynamics is ergodic, 
there is no need to worry about initial conditions, which might be said to be specified at a time t = — oo. 
Here, we are interested in how the equilibrium state is reached from a given initial state and must include 
into the action a term describing the initial preparation of the system. One has 

= S[0,0]+a[0,0] (3.2) 

where 0] p.ip describes the 'bulk' evolution of the system as derived from the Langevin equation 
while o"[0, 0] describes the initial conditions at time t = 0. As already pointed out by Mazenko ^7\, it 
may be written as 

a[(f),^] = --l drdr'0(O,r)a(r-r')0(O,r') (3.3) 
where it is implicitly admitted that (0(0, r)) = and a{r) is the initial two-point correlator 

a(r) := ^(0, 0; r + r', r') = (0(0, r + r')0(O, r')) (3.4) 

From spatial translation-invariance, it follows that a{r) = a{—r) which we shall admit throughout. 
We call the theory described by the action 5*0 alone the noiseless theory. 

For a free field, the noiseless and the thermal parts of the MSR action read (we also have set = 1) 

^o[0,0] := jdtdr^(^^-A<j) + v{t)<P 
5[0,0] := -T Jdtdr^^{t,r) 

S[0,0] = 5o[0,0]+5[0,0] (3.5) 

We shall refer to the contribution described by S and a as the thermal and initial noise, respectively. 

^In the systematic terminology of |17|. this should be rather called a response operator, because cj) will become an 
operator in a canonical quantization scheme of the action. The notion of a response field should have been reserved to 
the canonically conjugate variable of (f). Since we shall not use the operator formalism here, we shall simply, but sloppily, 
talk about as a response field. 
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We point out that field-theoretic studies of critical dynamics use a different initial term, namely [HHl 



^c[0, 0] = ? / dr (0(0, r) - m{r)f (3.6) 

This specifies an initial macroscopic state with spatially varying order parameter (0(0, r)) = m(r) and 
spatial correlations decaying on a finite scale proportional to tq^^. Ageing at criticahty was studied in 
the 0(n)-model using the e-expansion with the initial term ctc jHHl 1121 CHI lEl CS] • However, the use of 
(Tc instead of a for temperatures below criticality would lead to contradictions. 

Treating the noisy Langevin equation (jl.8|) as the classical equation of motion of the field-theory 
(MSR) action eqs. ()3.2I3.5I3.3|) has the following advantages. 

• Thermal fluctuations and initial conditions are explicitly included. To emphasize this, notice that 
the noisy contributions to the MSR action are 

S[(j),(f)] — So[(j),(j)] = — dudr du dr' (f){u,r)K{u,u';r — r')(f){u' ,r') 



= T6{u-u')6{r) + ^6{u)6{u')a{r) (3.7) 

where k includes both the effects of thermal and of initial-state fluctuations. 

The response field describes the thermal noise as can be seen from the equations of motion 
derived from the free-field MSR action ()3.5|1 

= A0(t,r)-t;(t)0(t,r) + 2T0(t,r) 

90(t,r) 



dt 



A0(t,r)-t;(t)0(t,r) (3.8) 



The first equation ()3.8|1 reduces to the Langevin equation ()1.8|1 provided one makes the formal 
identification 

r/(t,r) = 2r0(t,r) (3.9) 

Therefore, at the classical level, the stochastic Langevin equation eq. ()1.8p is described by two 
deterministic equations. In our case they are both of Schrodinger type and with opposite masses 
for the field and the response field 0. 

Averages of any ra-point function built from the fields 0, can be expressed in terms of the 
functional integral 

F{0(t„r,),0(t„r,)}) = J P0P0 exp (-5[0,0]) F{0(t„ r,), 0(t„ r,)} (3.10) 
with the normalization (1) = 1. 

For example adding a magnetic perturbation 5Hmag = — hip io the Hamiltonian (jl.llj) and 
then computing the mean of the order-parameter to first order in h, the relation i?(t, s; r, r') = 
{(pit] r)(f){s, r')) of eq. ()2.8|1 is easily reproduced. 

The use of Martin-Siggia-Rose formalism makes the machinery of the field-theoretic renormalization- 
group available |E3l [T7| IHII but we shall not pursue this here. 
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It will be useful to split the calculation of averages into two steps as follows 

^F{(/.(t„r,),0(t,,r,)}) = ^F{0(t„r,),0(t„r,)}exp 0] - (t[0, 0]) )^ 

(F{0(t„r,),0(t,,r,)})^ = J exp (-S'o[0,0]) r,), 0(t,, r,)} (3.11) 

where the notation ()o (and more generally the index 0) refers from now on to averages of the non- 
fiuctuating theory. This allows to make use of the Schrodinger invariance of the noiseless theory. 



3.2 Analytical results for free fields 

In this section, we find both response and correlation functions for free-field Martin-Siggia-Rose theory, 
as given by eqs. ()3.2I3.3I3.5|) . 



3.2.1 Two-point functions without noise 

The free Martin-Siggia-Rose action S[(f),(j)] has a Gaussian structure. We shall write it as 

S[(p,^] = jdudrdu'dr' <!?{u;rfQ{u,u';r,r')^{u';r') (3.12) 

where $ = f ~ ) is the two-component field built from and 0, and $^ stands for its transpose. The 

\(P J 
kernel Q reads 



Q{u,u';r,r') = ^ 



6{u - u')5{r - r') (-A - d^) + ^v{u) 

6{u - u')6{r - r') (- A + du) + ^v{u) -2k{u, u'; r - r') 

(3.13) 

This peculiar form of the Lagrangian density is quite suggestive as regards the Galilei-invariance. When 
K = 0, that is in absence of noise, the quadratic form Q is antidiagonal. This is one way of presenting 
the Bargman superselection rules and leads in particular to (00) = (00) = which is a manifestation 
of Galilei-invariance of the noiseless system. The presence of noise just breaks this symmetry. 

In order to study systematically the role of the noise, we shall expand around the non-fiuctuating 
theory. The correlation functions Cq, Cq and the linear response function Rq are 

Co{t,s;r,r') = 
Co{t,s;r,r') = 

Roit,s;r,r') = ^©(^ " ^) (M^ - 3))-'/^xp (-1^^=^^) (3.14) 

This follows since the quadratic form Q is antidiagonal the field can only be coupled to and the 
result is just the bare propagator of the theory. It is clear that only the antidiagonality of Q is important 
to derive this result and the explicit free-field form ()3.12|) is not required for ()3.14j) to hold. 

These results fully agree with the Schrodinger-invariance prediction eq. ()2.12I2.11|) with the identi- 
fications X = x = d/2, M. = 1/2 and ro = (47r)~'^/^. 
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A further manifestation of the Gahlei invariance of the noiseless theory is the fact that 

(0^0^)0 = (3.15) 

n m 

unless n = m as is easily checked. This is a further example of the Bargman superselection rule (j2.13j) 
and will be important in what follows. 

3.2.2 Two-point functions in presence of noise 

We now find the same two-point functions in the presence of noise. 

We begin with the response functions R{t, s; r; r') which is found from ()3.14|) by averaging with the 
noiseless weight exp (—So[(f),( 



R{t,s;r,r') = r)0(s, r') exp (^j dudRdu' dR' (l){u, R)k{u,u'; R- R')(l){u',r')^'^ (3.16) 

where k is given by eq. ()3.7p . Formally expanding the exponential and taking the Bargman superselec- 
tion rule ()3.15|1 into account, only the term of lowest order remains. We thus find 

Rit, s; r, r') = R^it, s; r, r') = ^0(t _ ,) (4vr(t - s))-'^' exp (^-^t-ll^j (3.17) 

and we see that R is independent of the noise. 

Next, the order-parameter correlation function reads 

C{t,s]r,r') = ^0(t,r)0(s,r')exp (^j dudRdu' dR' ^{u, R)k{u,u'; R - R')^{u' ,r')^^ (3.18) 
Again expanding the exponential and using eq. ()3.15|) . a single term remains and we readily find 

C{t,s;r,r') = J dRdudR' du' k{u,u'; R- R')R'^\t, s,u,u';r,r', R, R') 
R^o\t,s,u,u';r,r',R,R') := (^(l){t,r)^{s,r')^{u, R)^{u , R')^ ^ (3.19) 

where R'^^ is a noiseless four-point function. 

Repeating the same arguments as before, it is an easy task to compute the two-point correlations of 
response fields. They are 

C(t,s;r,r) = (3.20) 

We emphasize that eqs. ()3. 1613.1913.201) will hold for any theory satisfying the Bargman's superselection 
rule ()3.15j) . We shall come back to this in section 4. 

The four-point function R^^^ is simple to access for free fields since it factorizes into a product of 
two-point functions because of Wick's theorem. We have 

R'^\t, s, u, u'; r, r', R, R') = Ro{t, u; r, i^)i^o(s, u'; r', R') + i?o(t, u'; r, R')Ro{s, u; r', R) (3.21) 
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Together with the exphcit form of k, this yields the final result 

C(t,s;r,r') = Cth{t, s;r,r') + Cpr{t, s;r,r') 
Cth{t,s]r,r') = 2T jdudy Ro{t,u;r,y)Ro{s,u]r',y) 

Cprit,s;r,r') = J dy dy' Ro{t,0;r,y)a{y - y')Ro{s,0;r' ,y') (3.22) 

where we separated C in a thermal term Cth and an initial term Cpr- We clearly see that while the only 
contributions to C come from the noise, R does not depend on it. 

We summarize the results obtained so far as follows: 

1. It is satisfying that the well-known result C = {(fxp) = is naturally reproduced, see ()3.20j) . 

2. The independence eq. ()3.17|) of the two-time response function of T and of the initial correlations 
goes beyond the usual scaling arguments as reviewed in section 1. This explains to some extent 
the success of the existing confirmations of that prediction of local scale invariance. 

3. We arrive at an explicit expression for the two-time correlators, which are obtained in terms of 
a contraction of two response functions. We also see that the earlier result C = comes from 
neglecting both initial and thermal fluctuations. 

It is useful to present these results also in momentum space. Using spatial translation invariance, 
define the Fourier transform of any two-point function A{t, s; x, y) = A{t, s;x — y) as 

A{t,s]q):= / dr A{t,s;r)exp{-iq ■ r) (3.23) 




R{t,s;q) = ^exp{-q\t-s)) e{t-s) 

C,,{t,s;q) = 2T fdu^-^^exv{-q\t + s-2u)) 
Jo k\u) 

^ ^ ^ k(t)k(s) 

Cpr{t,s;q) = Ro{t, 0; q)a{q)Ro{s, 0; -q) = a{q) ^^^^^ ' exp {-q\t + s)) (3.24) 

where in the second line, the convention s < t has been used. In particular, if f (t) = the two-point 
correlation function takes an especially simple form C°(t, s; q) where 



C\t,s-q) 



T 

Hq) - — 



exp {-q\t + + ^ exp {-q\t - s)) (3.25) 



While both response and correlation functions depend on both t and s and therefore describe an ageing 
behaviour, there is an equilibrium regime 1 <^ t — s <^ s,t where we have the following simple expressions 

Req{t, s; q) = exp {-q'^{t - s)) 
a°q(t, 9) = ^ exp {-q'{t - s)) (3.26) 
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More generally, it is not difficult to show that Ceq(t, s; q) = {T/q^) e-'j'(*-^) (1 + ((sq2)~i, (t - s)/s)). 
In any case, we recover the fluctuation-dissipation theorem TReq^it, s; q) = dCeq{t, s; q)/ds in the equi- 
librium regime as it should be. 

Motivated from studies in spin glasses, it is sometimes attempted to separate correlation and response 
functions into an equilibrium and an 'ageing' part, viz. C = Ceq + Cage, R = Req + -Rage- In our case, 
we would have 

^age(t, s- q) = - 1^ exp {-q\t - s)) 

C°ge(t, q) = {a{q) - ^ exp {~q\t + s)) (3.27) 

For v{t) = 0, there is no 'ageing' part in the response function. 

From these expressions, we can already extract a few general properties of the ageing process. First, 
for systems quenched to below their critical temperature Tc > and described by a MSR Gaussian 
action ()3.1|) . it is known from the dynamical renormalization group that the flnal temperature T < Tc 
is an irrelevant parameter and furthermore T — > under renormalization j^l El • Then the long-time 
dynamics should be driven by the initial fluctuations, in agreement with eq. (|3.27|) with T = 0. On 
the other hand, for a critical quench T = T^, the situation is different in that both initial and thermal 
fluctuation may contribute to the long-time dynamics. From eq. ()3.27p we expect that the small-q 
behaviour of the term a{q) — Tc/q^ will determine the long-time dynamics. 



4 Consequences of local scale- invar iance of noiseless theories 
4.1 MSR formulation 

We now generalize the results of the previous section to any theory whose noiseless MSR action is 
Schrodinger-invariant. As we shall stay within the context of classical fleld theory, the symmetries of 
the MSR action are the same as for the corresponding Langevin equation (jl.8p . The formulation of 
the MSR action, using a non-conserved order parameter described by model A dynamics is almost 
unchanged with respect to section 3. We have 

0] = 5o[0, 4>] + S[<i), 4>] + ^[0, 0] (4.1) 

where 0] and cr[0, 0] are given by eqs. (j3.3l3.5|) and 

S„[^.?l=/d.d<?(f + ^) (4.2) 

We shall assume throughout that the noiseless action Sq is Schrodinger-invariant and that it includes 
an external time-dependent potential v = v{t). For a local effective potential Ti, it is known that spatial 
translation-invariance, dilatation-invariance (or dynamical scaling) and Galilei-invariance are sufficient 
for having Schrodinger invariance Therefore, from now on the dynamical exponent z = 2. 

In order to study the effects of the noise, we flrst discuss what becomes of the Galilei invariance of 
the noiseless theory. If the order parameter and the response fleld are quasiprimary, they should 
transform under a Galilei transformation t \—>- 1' = t and r r' = r — vt as (see eq. ()2.5|) ) 

0'(t',r') = Ut,r)4>{t,r) 

0'(t',r') = f-\t,r)<P{t,r) (4.3) 



13 



where the companion function reads [70] 

r) = exp 

The noisy contributions to the action transform as 



Mr ■ V — —v'^t 
2 



(4.4) 



dr dr'a{r - r')0(O, r)0(O, r') (/"^(O, r)f-\0, r') - l) (4.5) 



Therefore, for fixed temperature and initial conditions, the noise always destroys Galilei invariance.^ 
Thus, the only dynamic symmetries of the noisy Langevin equation should be space translations and 
dilatations (and possibly space rotations). 

Consequently, one would merely expect the following scaling forms 

t {r-r'f^ 
s ' if: — s 



C(t,s,r,rO 
R{t,s-r,r') 



(4.6) 



where Qc,r are undetermined scaling functions. For free-field theories, one would have a = h = d/2 — 1 
by dimensional counting. 



4.2 Response and correlation functions for the fluctuating theory 

The form of the scaling functions Qr^ will now be determined through a generalization of the expansions 
carried out in section 3. The main results will be eqs. ()4.8p and ()4.12|) . ()4.16|) . respectively. 

4.2.1 The two-point response function 

Consider a system which initially is at thermal equilibrium with temperature Tj (which fixes the initial 
correlator a{r)) and quench it at time t = to the final temperature T = Tf. The response function 
is still given by eq. ()3.16|) and the expansion of the exponential goes through as before. Because of the 
Bargman superselection rule (j3.15p which holds because of the assumed Galilei invariance of 5*0 only 
the lowest term survives and we obtain 

R{t,s;r,r')=Ro(t,s;r,r') (4.7) 

The expression of Rq has been derived earlier and using the gauge transform ()2.4j) we recover exactly 
the same result as in eq. ()2.12|) . namely 

R it, s; r, r') = 6.r. r, Q{t - s) {t - ^ exp [-^^I-llT] (4.8) 

k[s) V 2 [t-s) J 

We stress that, given only Galilei- and scale-invariance of the noiseless theory, this result should hold 
for any initial and final temperature Tj and Tf. At this stage, nothing has yet been said on the time- 
dependence of v{t). 

^The breaking of Galilei-invariance through thermal noise can be visuahzed as follows: consider a system in contact 
with a thermal bath at constant and uniform temperature T > 0. If the system moves with respect to the bath with a 
constant speed v, the apparent temperature measured in the system will depend on the angle between the direction of 
measurement and v. 
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4.2.2 The two-point correlation function 

The two-point correlation function is found from eq. ()3.19|1 . Again, the arguments of section 3 go 
through and we find, using the explicit expression ()3.7p for the kernel n 

C{t,s;r,r') =T J dudRR^Q\t,s,u;r,r',R) + ^ j dRdR' a{R- R')R^Q\t, s, 0,0; r,r', R, R') (4.9) 

(3) 

where Rq is the following three-point response function 

R^^\t,s,u;r,r',R) := (^<j){t;r)(f){s;r)^{u; R)^) ^ (4.10) 

and i^Q^"* was already defined in eq. ()3.19|) . This central result will be the basis for all what follows. 
Consequently^ the calculation of C requires the computation of the noiseless three- and four-point 
functions R^ and -Rq^^ This cannot entirely be done, since a general expression for R^^ is not yet 
available. Of course, one might hope that through an extension of the Schrodinger algebra sci)^ to some 
infinite-dimensional Lie algebra techniques analogous to 2D conformal invariance |3j might become 
applicable but the formulation of just such an extension is an open problem. 

Here, we shall restrict to the case of vanishing initial correlations, that is 

a{R) = aoS{R) (4.11) 

which corresponds to an infinite initial temperature Tj = oo and where is a normalization constant. In 
fact, from renormalization group arguments the long-time behaviour of any system which is prepared in 
the high-temperature or paramagnetic phase should be described by this initial condition [HI E] • Then 

C{t,s) = T jdudRR^o\t,s,u;R) + ^ jdRRP{t,s,0;R) (4.12) 
R^^\t,s,u;r) := R^i\t, s,u;y,y,r + y) = U{t;y)<P{s;y)4>{u;r + yf) (4.13) 





Here, the field 0^ is a composite field with mass —2M. and scaling dimension 2x2. Only for free fields, 
one has X2 = x. 

Now, the three-point function of a Schrodinger- invariant theory with v{t) =0 is well-known since a 
long time [JSlEn]- Denoting by TZ the response function in the case where v{t) = 0, we have 



T^Q^"* (t, s, u; r) = Tl''^^ (t, s, u) exp 



M t + s-2u , 



^ ( ^—^ 

{t — u){s — u) 



2 {s-u){t-u) 

n'i\t,s,u) = &{t - u)&{s - u) {t - u)-^' {s - u)-^' {t - sy^'^^' (4.14) 

where \1' is an arbitrary scaling function and is the 'mass' of the field (p. This result is brought to 
the case at hand through the gauge transformation ()2.4|) and we find 

R'o\t.s,u;r) = ^^^n^^\t,s,u;r) (4.15) 

Combining ()4.12|) and ()4.15p we obtain, for the first time, a generic prediction for the form of the two- 
point correlation function. It is remarkable that under the condition 1)4.1111 only the noiseless three-point 
response functions are required in order to predict any two-time autocorrelator. 
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It is useful to write down the autocorrelation function in the form C(t,s) = Cthit,s) + Cpr{t,s). 
Here the thermal part Cth and the preparation part Cpr are given by 



Cth{t,s) 



rpgd/2 + l-X-X2 



- 1 



X2—x—d/2 pi 



de 



k{t)k{s) 



pr 



2 P(0) 
dR exp 



d/2-x 



X2^x—d/2 



-^^ Wi- 



d/2-X2 



t/s+1- 29 
t/s-1 



- 1 



t/s + 1 
t/s-1 



Mw 



(4.16) 



and we have explicitly used s < t. 

As they stand, the above expressions for C{t, s) do not yet necessarily describe a dynamical scaling 
behaviour, since the form of k{t) is still completely general. We now assume in addition that we are 
dealing with a system with dynamical scaling. In order to reproduce the usual phenomenology of ageing 
systems, we must have, at least for sufficiently large times (see section 2) 



k{t) ~ kot^ 



(4.17) 



Comparing the general form of R{t, s) as given in eq. ()4.8|) with the phenomenologically expected scaling 
eqs. (jl.2ll.3p . we read off 

(4.18) 



R 



X = 1 + a , r = 1 + a — 

2 



and in particular, the scaling function p.7|) is recovered. On the other hand, for C{t,s) we find the 
following scaling form, written down separately for the thermal and the initial term, where y = t/s > 1 
is the scaling variable 



G,(t,.) = Ts-'^-fS'iy) 
Cpr{t,s) = s-'^^fl\y) 



(4.19) 



where the scaling functions and are given by 



y' (y-i) 



X2—x~d/2 



doe-'F [^y-e){i 



^d/2~X2 



y+l-2e 



^„d/2-S2+r (^y _ ^Y2-x-d/2 ^ ( 1/+ 1 



(4.20) 



with the non-equilibrium exponents 



Hh 



X + X2 - 1 - d/2 
x + X2-2r - d/2 



(4.21) 



Provided that $(1) is finite, the asymptotic behaviour of the scaling functions for y large can be worked 
out. We expect fciy) ~ y~^c:/'^ and find 



X'^ = A^J = 2(x - r) 
Therefore, comparing this with ()4.18|1 . we have shown: 



(4.22) 
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For any system with an infinite-temperature initial state l{4.11\ ), quenched to a temperature T < 
and whose noiseless part is locally scale-invariant with z = 2, one has 

Ac = A^, (4.23) 

For non-equilibrium critical dynamics (that is T = T^) the same conclusion can be drawn if after 
renormalization one still has z = 2. While eq. (|4.23p certainly agrees with the evidence available from 
models studied either analytically or numerically, we are not aware of any other general proof of this 
equality between the autocorrelation and autoresponse exponents for a fully disordered initial state. 

In sections 5-7, we shall present extensive tests of the prediction ()4. 1914. 2014. 21|) for C and ()4.8l4.18p 
for i?, respectively. The main hypothesis going into it is the requirement of Galilei-invariance of the 
noiseless theory, while the other conditions appear to be habitually admitted in the description of 
physical ageing. 



4.2.3 Autocorrelation function in phase-ordering kinetics 

In order to understand the result ()4.19|) for C(t, s) better, we now study the two contributions separately. 
First, we consider the 'preparation' part Cp^. This term is expected to describe the late-time behaviour 
of a system quenched to a temperature T < Tc. Indeed, renormalization-group arguments show jH] that 
in this case T is an irrelevant variable and is renormalized towards zero. Then the thermal contribution 
Cth vanishes. Therefore, the non-equilibrium exponent b = bpr is read off from eq. ()4.2H) . In addition, 
we know that 6 = in the low-temperature phase. This implies 

X2-x = ^-Xc<0 (4.24) 

because of a well-known inequality jHS]- Only for free fields, one has Ac = d/2, otherwise X2 is a new 
nontrivial exponent. In the scaling limit, we thus have C{t, s) = fc{t/s) where 

fc{y) = f y^-/^ {y - 1)-'^ $ (^) (4.25) 

The form of fc{y) still depends on the unknown function $(w) which in turn depends on \l/(p), see 
(|4.16p . We attempt to fix its form and reconsider the noiseless response function Rl^\t,s,0;r) which 
describes a response of the autocorrelation C{t, s) = (0(t)0(s)). It appears to be a reasonable require- 
ment that there should be no singularity in Rq when t = s. Using the explicit form eqs. ()4.14f4.15j) 
this leads to the following limit behaviour 

^(p) ~ ^oP^^"'^^^ ; P (4.26) 

where \l/o is a constant. If this leading term should be still accurate for larger values of p, the following 
expression for the scaling function ^{w) is found 

where 5*^ is the surface of the unit sphere in d dimensions. Eq. ()4.27p should hold in the w ^ oo limit. 
Provided this from is still valid for all w, we would obtain the following simplified form, with y = t/s 

0,,,.-^£(fe±i)!)-^"^Mi i^f'" ; ,4.28) 

where we also related the leading constant to the squared equilibrium magnetization M^^, in order to 
recover the usual scaling form C{t, s) = M^^fcit/ s) with /c(l) = 1, see [H^ . 
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4.2.4 Autocorrelation function for critical dynamics 



Second, let us turn to the thermal part. It should dominate the autocorrelation for quenches to the 
critical point T = T^, given the initial condition ()4.11|) and under the assumption that z = 2 even at 
criticality. From eq. ()4.21|) . the exponent h = hth can be read off and we have 

X2 = h-a + ^ (4.29) 

Under the stated conditions, the preparation term drops out at large times and we find 

C{t,s) = T,s-'fc{t/s) 

fc{y) = {y - if''"-' de r [{y-e){i- e)r' $ (^ ^ + (4.30) 

At criticality, one expects the following relationship between the nonequilibrium exponents a and h 

a = b='l = '-I±^ (4.31) 

uz z 

where (5,^,7] are well-known equilibrium critical exponents (see section 1). To understand eq. ()4.31|1 . 
recall that C{s,s) ~ s~^. On the other hand, from the space-time scaling \r\^ ~ t, one expects the 
equilibrium correlator to decay as Ceq ~ \r\~^^ and the second equality in eq. ()4.3H) follows. Finally, 
a = 6 is a necessary condition for having a non-vanishing limit fiuctuation-dissipation ratio X^o- 

Now, if we let a = & in ()4.3U|) . we find 

fciy) = y'{y- i)"^"" (^7^) ^^'^^^ 







which is the most general form compatible with the standard phenomenological constraints. An ap- 
proximate form of the scaling function $(w) may be obtained from the requirement that the three- 
point response function TZq {t , s , u; r) is non-singular for t = s. We now have X2 = d/2 and find 
\E'(p) ~ \E'o,c p'^"''^^ as p 0. This leads to, for w — * oo 

$(w) ^ $o,cW-^-" (4.33) 

where $o,c is a constant. If in addition we may use this form also for finite values of w, we would obtain 
the simplified form 

fciy) ~ %,c yi+'^-^^Z^ d^ QXc-2-2a + 1 _ 29)-'-'' ; T = T, (4.34) 

Jo 

Summarizing, the phenomenological comparison of the autocorrelation function, as predicted by 
Schrodinger-invariance, and assuming a totally disordered initial state, with simulational or experimen- 
tal data will be based on eqs. ()4.25|) and ()4.32|) for quenches to T < and T = T^, respectively. In full 
generality this will allow to obtain information on the scaling function $(w). If in addition the heuristic 

(3) 

idea of the absence of singularities at t = s in the three-point response function Rq should be sufficient 
to fix the form of this response function, the simplified forms ()4.28p4.34j) apply and the scaling function 
fc{y) is completely specified in terms of the exponents a and Ac- 
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5 Tests of local scale- invar iance in exactly solvable models 



In this and the next two sections we describe phenomenological tests of the predictions of Schrodinger- 
invariance, that is local scale- invariance with 2; = 2, which were derived in section 4 in concrete physical 
models of ageing behaviour. 



5.1 Kinetic spherical model 

The kinetic spherical model is often formulated in a field-theoretic fashion as the n 00 limit of the 
0(n)-symmetric vector model. For our purposes, it is more useful to start directly from a lattice system 
and to take the continuum limit later, following [^1^1^17711751 173]. 

Consider a hypercubic lattice with J\f sites. At each site r there is a real time- dependent variable 
(j){t, r) such that the mean spherical constraint 



(5.1) 



holds. The Hamiltonian is 7i = — X]<r r'> *^(^' ''")'^('''' ''"') "where the sum extends over pairs of nearest 
neighbour sites. The (non-conserved) dynamics is given in terms of a Langevin equation of the type 
eq. (HH) 

= 5^0(t,n)-(2rf + .;(t))0(t,r) + r/(t,r) (5.2) 

n{r) 

~ A(p{t,r) - v{t)(t){t,r) +r]{t,r) (5.3) 

where n(r) runs over the nearest neighbours of the site r. In the second line, a formal continuum limit 
was taken (for simplicity, all rescalings with powers of the lattice constant a were suppressed). Finally, 
ri{t,r) is the usual Gaussian noise (see section 1) and the Lagrange multiplier v{t) is fixed such that 
the mean spherical constraint is satisfied.^ Therefore, the kinetic spherical model perfectly fits into the 
context of a Schrodinger equation in a time- dependent potential v = v{t) discussed in section 3. We 
therefore expect that the free-field predictions eqs. ()3.17I3.22|) will be fully confirmed. 

In order to see this explicitly, we recall the elements of the exact solution of the Langevin equation 
(Q. If we set 

^(t)=exp(2 / duv{u)] (5.4) 







it can be shown |22linHllZZ! that g{t) is the unique solution of the Volterra integral equation 

g{t) = A{t) + 2T fdt' fit - t')g{t') (5.5) 
^0 



where g{0) = 1 and 



fit) = e(t) (e-%(4t))' 



1 



A{t) = I dq yXr)e-'-''^''^''- (5.6) 



careful study shows that provided the Umit A/" ^ cxd is taken before any long-time limit, the mean spherical 
constraint (|5.1|l and a full, non-averaged, spherical constraint lead to the same results (29) . 
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and uj{q) = 2 — cosgj) is the lattice dispersion relation, B the Brillouin zone while Jq is a modified 

Bessel function. We stress that because of eq. fl5.5|) . g{t) does not depend on the order parameter (f){t, r). 

We now show that the prediction ()H.17|H.22j) for the two-point functions can be fully reproduced. 
We begin with the response function. In the spherical model, it is exactly given by jHHl EZl 



R{t,s;r,r') 



6{<P{t,r)) 



6h{s,r') 



h=0 



ne-^(*-^)/..-.; (2(t-.)) 



9it) 



(47r(t-s))-'^/2exp 



r\2 



{r - r') 

4(t - s) 



exp 



duv{u) e(t - s) 



(5.7) 



where in the second line the limit t — s^l was taken and Ir is a modified Bessel function. In particular, 
this reproduces the known autoresponse function 38^ 



Rit, s) = Kit, s; r, r) = /((t - s)/2)y/gis)/git) 



{5.i 



This is in exact agreement with eq. ()3.17|) and we identify the mass Ai = 1/2. 

We now turn to the correlator but the sake of brevity merely deal with the autocorrelator explicitly. 
In the spherical model one has (with t > s) jHHl EZl US] 



C(t,s) = (0(t,r)0(s,r)) 



1 



(5.9) 



In order to show how to recover this explicitly from eq. ()3.22|) . we write again C{t, s) =: Cpr + Cth and 
discuss the two terms separately. The first one is, where we use the explicit form (j5.7|) oi R = Rq 



pr 



dydy' (4vrt47rs)-'^/%/4^4^exp 



9{t) g{s) 



[r - yf {r - y'f 



At 



As 



(2vr) 



-d/2 



f dqdq' e-'?^*-'?'^-+i(9+9')- . 



where 



dydy' a{y - y') e 



(5.10) 



(5.11) 



and we used g{Q) = 1. In order to calculate J^, we set C = y — y', C = y + y'- Then the Jacobian 
\d{y,y')/d{C,C)\ = 2"^^ and we have 



J, = 2-^ / dCdC'a(C)e 

= {2TxY5{q + q') [ dCa{C)e-"^< 



We finally obtain 



a 



pr 



\-. [ dqe-'''(*+^) / dCa(C)e~''^-^ 



(5.12) 



(5.13) 
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When the support of a(C) is restricted to the hypercubic lattice, we therefore have indeed for long times 
t + s > 1 

C,. . ^^^^±m (5.14) 
in agreement with the first term in ()5.9|) . The second term in ()3.22|) is analyzed as follows, for t > s 



th 



poo 

2T / du 
Jo 



X exp 



dy (47r(t — u) 4it{s 
(r 



u)) 



-d/2 9{u)g{u) 



y 



y) 



A{t - u) 4(s - u) 



0(t-M)0(s 



2T / du 
Jo 

2T 

V9it)9is) 

2T 
V9it)9is) 



dqdq'dy e-'j'*-9'^*+{q''+q")"+i(q+q') »' g-i(9+9')-y 



9{u) 



dug{u) (27r 



dq e 



-2q'^{it+s)/2-u) 



dug{u)f 



t + s 



— u 



(5.15) 



where the last line holds for sufficiently large arguments of the function /. Taken together with Cp^, 
the expected agreement between the general result eq. ()3.22|) and the exact expression ()5.9|) for the 
spherical model is thus recovered. The full space-time correlator can be checked similarly. 

It is interesting to note that for the confirmation of i?, see fl5.7|) . we need the condition t — s ^ 1, 
while for the confirmation of C, we also need s ^ 1 (which implies t + s ^ 1). 

So far, we have not yet used the explicit form of gii) which follows from eq. ()5.5p . Indeed, it is 
well-known that for long times, one has I77| HHl 

git) ~ (5.16) 

For the fully disordered initial conditions ()4.11|) . the exponent F takes the following values: (i) for 
T < Tc, one has F = d/A and (ii) for T = T^, one has F = 1 - if 2 < d < 4 and T = if > 4. 
This is exactly the form expected from a potential of the form t>(if:) ~ F/t, see eq. fl2.14j) . 

It is instructive to compare also the explicit result for the two-point functions with the expectations 
coming from local scale-invariance. For the two-time response function, this confirmation has already 
been carried out for both the autoresponse function R{t, s) as well as the space-time response R{t, s; r) 
|45| Hni EZl and need not be repeated here. Therefore we concentrate on the two-time autocorrelation 
function C{t, s). First, we consider the case T < Tc where from the exact solution it is well-known that 



Cit,s)=Mlfcit/s) , fc{y) 



Ay 



-d/A 



(5.17) 



and we read off from the y — > oo limit the exponent Ac = d/2. Clearly, this exact result is in full 
agreement with the prediction ()4.28|) of local scale invariance. Second, we consider the case T = T^. 
Then the exact solution gives 38 

,,,,, , ^,{y - i)-'/'^V-'/\y + 1)-' ;if2<d<4 

_2_[(^_l)-rf/2+l_(^^l)-d/2+lj .ifd>4 

{5.1i 
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From this, we read off the exponents a = h = d/2 — 1 for dX\ d > 2 and furthermore Ac = 3c//2 — 2 
for 2 < d < A and Xq = d for d > 4, respectively. Inserting the exponent values into ()4.34|) it is 
straightforward to show that eq. (jS.lSj) is indeed reproduced. 

In particular, we see that in the spherical model the asymptotic forms eqs. ()4.27l4.33p . respectively, 
are indeed exact as should be expected for a free-field theory. 

In conclusion, the exact results of the kinetic spherical model for both two-time correlation and 
response functions are in full agreement with the predictions of Schrodinger invariance. 



5.2 XY model in spin- wave approximation 

Another system which can be exactly analysed is the kinetic XY model in spin-wave approximation. 
As we shall see, it provides an instructive example on the correct identification of the quasi-primary 
scaling fields in a given model. 



5.2.1 Formulation and observables 

The XY model describes the interaction between planar spin variables 

5(r) = cos(0(r))ei + sin(0(r))e2 = ( ""^tiy) ) ^^'^^^ 

which are attached to the sites r of a c?-dimensional hypercubic lattice and </'(r) is the phase. The 
Hamiltonian is 

n\<P] = -J2 Sir) ■ Sir') = - 5^ cos (0(r) - 0(r')) (5.20) 

where the coupling constant J has been set to unity and the sum runs over nearest neighbours. The 
relaxational dynamics is assumed to be described by a Langevin equation. We prepare the system 
initially a temperature Tj and quench it at time t = to the final temperature T = Tf, so that the 
angular variable obeys [0] 

^-^-'-)- 

where t] represents an uncorrelated Gaussian noise with zero mean and variance 

{r]it, r)r/(t', r')) = 2Tj5(t - t')6ir - r') (5.22) 

Here, we shall exclusively study the coarsening dynamics in the low-temperature regime, that is 

T,,Tf<^T,id) (5.23) 

where Tdd) is the critical temperature of the XY model in d dimensions (if d = 2, Tc(2) = Tkt is the 
Kosterlitz-Thouless temperature of the transition). Then the so-called spin- wave approximation [S^ll] 
can be used which amounts to expand H in powers of 0(r) — 0(r'). Shifting the energy by a constant, 
the Hamiltonian reads to lowest order 

n[<P] = 1 fdr (V0(r))2 (5.24) 
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In writing this, we have imphcitly absorbed the spin-wave stiffness jHU El] into a redefinition of the 
temperatures. In 2D, it is known that below Tkt, any vortices present will be tightly bound and for 
distances larger than the characteristic pair size, the XY model renormalizes to the Hamiltonian ()5.24p 



We are interested in the properties of the two-point functions. It appears natural to define two-time 
correlation and linear response functions in terms of the magnetic variables 

V{t,s-r,r') := r) ■ r')) = (cos r) - r'))) 

p(t,.;r,r') := hm (5.25) 
h-^Q dh{s,r') 

where the response is found by adding a term 6Htnag = J2r^ ' ^'^ Alternatively, one may also 
consider the analogous quantities defined for the angular variables 

C{t,s-r,r') := (0(t, r)0(s, r')) 

R{t,s;r,r') := Jim (5.26) 

h-^o dh{s, r') 

(where a perturbation (57iang = Ylir ^'^ should have been added). 

In order that the spin-wave approximation be applicable, we must start from an (almost) ordered 
initial state of the system. Therefore, we require the following initial value for the magnetic correlator, 
which reads in Fourier space 

a{q) = d{0,0;q) = ^-^ = -, (5.27) 

where fjiTi) is the standard equilibrium critical exponent and the the relation 27Tri{Ti) = Ti valid in the 
spin- wave approximation was used, see [Ell IHH Hj • 



5.2.2 Non-equilibrium statistical field theory 

As before, we introduce a Martin-Siggia-Rose formalism which characterizes the ^ystem is term of an 
action S[(f), (j), h] depending on the phase field and an associated response field (j) and we also include 
a (possibly space-dependent) magnetic field h = J2i hiCi. We decompose the action 



5[0,0] = S[0,0,/i]+a[ 



(5.28) 



into a bulk term 0, h,] and an initial term cr[0, 0]. These two terms include the thermal and the 
initial noise. Explicitly 



dt dr 



— A0 -|- sin(0)/ii — cos{(j))h2 



T dtdr (j){t,r)(f){t,r) 



(5.29) 



and 



i /drdr'0(O,r)a(r-r')0(O,r') 



(5.30) 



where the function a{r) describes the initial conditions according to eq. ()5.27|) . 

'''For quenches from above Tkt in 2D, vortex configurations also become important and this leads to logarithmic 
scaling, see [SH] for details. 



23 



We now simplify the general expressions for the two-point functions. There is nothing to do for the 
angular correlation function C{t, s; r, r') = {(f){t, r)(j){s, r')) and we start with the magnetic correlation 
function T which is given by 



T{t, s; r, r') = / V(f)V(j) cos r) - (f){s, r')) exp ( -S[(f) 



(5.31) 



For a vanishing magnetic field the bulk action S[(/), (f), 0] is a quadratic form in the fields 0, (f) which are 
therefore Gaussian. Standard techniques explained in appendix A lead to 



Tit, s] r, r') = exp 
= exp 



Cit,s;r,r')- 



C{t,t-r,r) + C{s,s;r',r') 



(5.32) 
(5.33) 



where in the second line the argument of the exponential was expanded. This gives F in terms of angular 
correlators C. 

Next, we consider the response functions. For the angular response R, we quote from the MSR 
formalism the standard result 

/?(t,s;r,r') = (0(t,r)0(s,r')) (5.34) 

It remains to consider the response of the spin vector S at time t and position r to some magnetic field 
h{s, r') at time s and position r'. From the definition ()5.25j) we have 



p{t, s] r, r') = lim 



((cos0(t,r))) - (cos0(t,r)) ((sin0(t, r))) - (sin0(t, r)) 



hi{s,r') 



hois, r' 



(5.35) 



where the average ((■)) is to be taken with a magnetic field. Expanding the action (j5.28|) to first order 
in both components hi, h2 of the magnetic field, we find 



((cos0(t, r))) = (cos0(t, r)) + hi{s, r')(cos0(t, r) sin0(s, r')(j){s, r')) 
((sin0(t, r))) = (sin0(t, r)) — h2{s, r')(sin0(t, r) cos0(s, r')(f){s, r')) 

It follows that the response function can be expressed as 

pit, s; r, r') = (0(s, r') sin (0(t, r) - (f){s, r'))) 

Since both fields 0, (p are Gaussian, it can further be shown that (see appendix A) 

p{t,s;r,r') = R{t,s;r,r')r{t,s;r,r') 



(5.36) 



(5.37) 



(5.38) 



and we see explicitly that the relationship between R and p is non-trivial. That no higher correlators 
than the magnetic two-point correlation function F enter is a consequence of the Gaussian nature of 
the theory at hand. 

Eqs. ()5. 3315. 3415. 38|) are the main results of this subsection. Before we can evaluate them, we need 
some information on the validity of the spin-wave approximation. 
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5.2.3 Remarks on the validity of the spin- wave approximation 



We need a criterion informing us up to wliat point tlie results on T and p derived in tlie previous 
subsection within the spin-wave approximation should be reliable. 

The correlation function C{t,s;r,r') has already been obtained above and is given in eq. (|3.25p . 
For our choice ()5.27|) of initial conditions, its Fourier transform C{t, s; q) is [2111 



C(t, s- q) = (T, - Tf)G{t + s-q)+ TfG{t - s; q) (5.39) 

where G is given by 

q) ■= \ exp {-q\u + A^)) (5.40) 

and we have explicitly introduced an uv-cutoff which simulates the lattice spacing (we shall let A ^ 
at the end). Therefore, a two-point correlation function (00) is of order 0(Tj,Tj) whereas a response 
function (00) is of order 0(1) in the initial and final temperatures. 

In order to discuss further the validity of the spin- wave approximation, we keep the next-order term 
as well and consider the Hamiltonian 

n<P] ='H, + \jdr [(V0(r))2 + g, (V0(r))^] (5.41) 

where is some constant. A straightforward calculation shows that to first order in 5(4, the correction 
to the spin-wave approximation of the two-point correlation function is given by 

5G{t,s-r,r') = y"dndi2(0(t,r)0(s,r')0(n,i2)V(V0(«,i?))')^o|^(^^^ ' (t^) ) ^^'^^^ 

where the six-point function is factorized into two-point function by Wick's theorem. As a result, the 
spin-wave approximation is a first-order approximation in the initial and final temperatures. Consistent 
expressions of two-point functions must be expanded to first order in Ti,Tf. Higher-order terms in Tjj 
calculated within the spin-wave approximation should not be expected to be reliable. 



5.2.4 Correlation and response functions in the spin-wave approximation 

Finally we are ready to list the result for two-time correlation and response functions in the spin-wave 
approximation. From the previous subsection we know that G must be of first order in Ti,Tf. The 
consistent result for the magnetic correlation function is therefore obtained by expanding eq. ()5.33|) to 
first order in temperature. Then 

T{t, s; r, r') ^ 1 + G{t, s; r, r') - ^ [Git, t; r, r) + C(s, s; r\ r')] (5.43) 

A more suggestive form of this is found as follows. We have 

(5(t, r)) ■ {S{s, r')) = (cos 0(t, r)) (cos 0(s, r')) + (sin 0(t, r)) (sin 0(s, r')) 

^ 1 - ^(0(t,r)2)(0(.,r')^) + (0(t,r))(0(.,r')) + ■ ■• (5.44) 

where in the second line we performed a low-temperature expansion which must be kept to second order 
in in order to be of first order in the temperature, since G = (00) = 0(T). Furthermore, because 
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of the 1-^ —0 inversion symmetry, (0) = 0. Inserting this into (j5.43j) we find, of course only in the 
context of the spin-wave approximation 

{Sit, r) ■ Sis, r')) - {Sit, r)) ■ r')) = Cit, s; r, r') (5.45) 

and the relation between F and C is finally clarified (the equilibrium version of this is well-known, see 
[Sg, sect. 4.2.2). 

For notational simplicity, we shall now concentrate on the autocorrelation and autoresponse func- 
tions. First, the angular correlation function is given by eqs. ()5.39|) and ()5.40j) . We have 



2(47r)-''/2 
d-2 



[(Ti-T/)(t + s + A2)i-'^/2 + T/(t-s + A2) 1-^/2] ; if > 2 



Cit,s)={ (5.46) 
(47r)-i [iTf - Ti) ln(t + s + A^) - T/ ln(t - s + A^)] ; if d = 2. 

Second, the autoresponse functions are given by 

pit,s) = Rit,s)(^ + Cit,s)-]^[Cit,t) + Cis,s)] 

Rit,s) = [AT^it - s + A^)]""'^ (5.47) 
These results require a detailed discussion. 

1. The 2D XY model was studied in detail in the spin- wave approximation before |3] and we 
now show that their results, although they might at first sight appear to be different, agree with 
eqs. fl5.43p5.46|5.47j) . For notational simplicity, we restrict to Tj = 0, = T. First, the magnetic 
autocorrelation function is (4j eq. (11)] 

Tit,s) = {Sit,r)-Sis,r)) ' ^ 



(2t + A2)(2s + A2)(t-s + A2)2 



exp 



viT)^^f A\t + s + A 



2\2 



(2t + A2)(2s + A2)(t-s + A2)s 



(5.48) 



1 + ^ (\nit + s + A^)- ln(t - s + A^) - 1 ln(2t + A') - ^ ln(2s + A^) + In A^ ) + 0(T'^ 



since the spin-wave approximation is only consistent to lowest order in T. It is now clear that the 
above result is reproduced by inserting C(t, s) from eq. ()5.46|) with d = 2 into ()5.43|) . Second, the 
linear spin response in 2D is jH eq. (13)] 

, 6{Sit,r)) 1 f A4(t + 5 + A2^2 ^viT)/4 

pit, s) = lim — ' 



h^o 6his,r) 47r(t-s + A2) V(2t + A2)(2s + A2)(t - s + A2)2^ 
in agreement with eq. (j5.47j) with ci = 2, as it should be. 

In 2D, the results for r(t, s) and pit, s) were confirmed by a simulational study with an ordered 
initial state and < |lj. 

2. For Ti = Tf, the two-point functions are stationary. This is only to be expected, since in this case 
the system was prepared in an equilibrium state and remains there. 
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3. Ageing occurs when Tj 7^ Tj. Time-translation invariance is broken and we proceed to analyse the 
resulting scaling behaviour. Now, there are in principle two equally appealing sets of variables. 
First, we may choose to work with the angular correlation function C and its associated response 
R. Recalling the scaling forms introduced in section 1 (see eqs. (I1.2|1.3|) ) we shall characterize 
them by the exponents a, b, Xc, Xr. Second, we may prefer instead to work with the magnetic 
correlation function F and its associated response p. We shall use the same scaling forms, but 
for clarity we shall denote the corresponding exponents by a', b', Xq, X'j^. These exponents are 
straightforwardly read off in the ageing regime where t, s,t — s ^ A and we collect the results in 
table [U The exponent = {rji + rjf)/2 was already known 

Table 1: Ageing exponents of the (i- dimensional XY model in the spin- wave approximation. Here 
rji J = ri{Ti f) = Tjj/(27r) describe the initial and final correlation exponents. 



angular correlation and response 




a b Xc Xr 


d = 2 
d>2 


2 
d/2-1 d/2-1 d d 


magnetic correlation and response 




a' b' X'c X'ji 


d = 2 
d>2 


r^f/2 r^f/2 (Vi + Vf)/^ 2 + (r/, + r/y)/2 
d/2-1 d/2-1 d d 



We see that the exponents satisfy the equalities a = b and a' = b' expected for nonequilibrium 
critical dynamics and point out that in 21?, the autocorrelation and autoresponse exponents are 
different: Xr — Xc = ^r — = 2. This effect comes from the non-disordered initial condition of 
the spins, as explained first in [77] . 

4. Having discussed the values of the ageing exponents, we wish to compare the form of the scaling 
functions with the predictions of local scale-invariance as derived in section 4. This requires, 
however, the correct identification of the quasiprimary fields in our system, see section 2. It is 
only the quasiprimary fields which are expected to transform in a simple way under a local scale- 
transformation and the transformation laws of more complicated fields built from quasiprimary 
fields must derived accordingly. 

In the case at hand, it is clear from the complicated structure of the magnetic correlation and 
response functions that the magnetic order parameter S'(t, r) does not correspond to a quasipri- 
mary field. Rather, the quasiprimary field should be identified with the phase (/)(t, r). Indeed, the 
form of the angular response s) is in perfect agreement with the prediction ()4.8l4.17j) of local 
scale-invariance. This suggests that the response field 0(s, r) should be quasiprimary as well. 

5. Having thus identified and as quasiprimary fields of the model, it is now clear that the angular 
autocorrelation function C(t, s) should be compared to the critical dynamics correlation function 
as derived in section 4. However, a direct comparison with eq. fl4.34|) is not possible, since in its 
derivation fully disordered initial conditions were assumed. 

We shall therefore proceed in two steps. First, we shall consider the case Tj = 0. Because of 
our initial conditions ()5.27|) the initial correlator then vanishes and eq. ()4.34|) should now hold 
true. Second, we shall show that in the context of the free-field theory underlying the spin-wave 
approximation of the XY model, the restriction to uncorrelated initial states can be lifted. 
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We now set Tj = and Tf = T. From the exponents in table HJ the predicted autocorrelator 
scahng function follows from ()4.34|) as 



fc{y) = <f fdO {y + 1- 29)-'^' ■ d>2 (5.50) 

and we see immediately that this is in agreement with the explicit angular correlator ()5.46jl . upon 
identification of $o- 

Finally, if we also allow for > 0, there is a contribution to C{t,s) from the initial condition. 
We then return to the basic result ()4.9|) and decompose C(t, s) = Cth + Cpr- The thermal term 
Cth was treated before and the preparation term is analysed in appendix B, with the result 

Cpr='^^^^^Ut + sf-'''^ ■ d>2 (5.51) 
in complete agreement with ()5.46p . The case c? = 2 is treated similarly. 



In conclusion, the two-time autocorrelation and autoresponse functions of the XY model treated 
in spin-wave approximation are in perfect agreement with local scale-invar iance, provided the angular 
variable and its associated response field are identified as the quasiprimary fields of the model. 



5.3 Fluctuation-dissipation relations in the XY model 

Having checked that both correlation and response functions agree with the local scale-invariance predic- 
tion, we now inquire what can be said on the approach of the model towards equilibrium. A convenient 
way to study this is through the so-called fluctuation-dissipation ratio [21], see section 1. Since we have 
seen that in the XY model angular and magnetic observables behave quite differently, it is convenient 
to define two distinct fluctuation-dissipation ratios, namely 

X{t,s) := T,R{t,s)(^-^^\ ' (5.52) 

Of particular interest will be the limit fluctuation-dissipation ratios defined in eq. p.5j) and similarly 
Soo. We have seen before that the scaling of autocorrelation and autoresponse functions is according 
to the expectations of nonequilibrium critical dynamics. In this case, according to the Godreche-Luck 
conjecture [H3, X^o and Soo should be universal. 

We shall use the available exact results in the XY model to test this conjecture by studying the 
dependence of X and S on the ratio a := Ti/Tf oi initial and final temperatures. 



5.3.1 Fluctuation-dissipation ratio for magnetic variables 



The fluctuation-dissipation ratio S(t, s) obtained from the magnetic correlation and response functions 
reads, with y = t/s 



1+1- 



Tf 



y-1 
y + l 



d/2 



y-1 



d/2 



(5.53) 
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For Tj = and d = 2 this was already known jl] . In the quasiequihbrium regime y — 1 the fluctuation- 
dissipation theorem should hold. Indeed we find lim^^^i = 1 which confirms that S should be 
well-defined. For large values of y, that is for well-separated times, we have 



and therefore Soo = 0, indeed a universal constant. It is remarkable that the asymptotic value of 
should be independent of d and that is agrees with the value Sqo = of phase-ordering kinetics in 
the low-temperature phase with an ordered, non-critical equilibrium state. This kind of result should 
be more typical of an ordered ferromagnetic equilibrium state as it occurs for d > 2 but is not really 
expected for d = 2 since the equilibrium 2D XY model is critical even below Tkt- 

Finally, for large y the asymptotic form of is independent on whether the system is cooled or 
heated. The temperatures merely enter into a scaling amplitude. 

5.3.2 Fluctuation-dissipation ratio for angular variables 

In the same way, the fiuctuation-dissipation ratio for the angular variables is found. It reads 



As before, in the quasiequihbrium regime t ~ s, s) = 1. Surprisingly, however, for large values of 
y = t/s, the limit fiuctuation-dissipation ratio 



depends continuously on a = Ti/Tf. We recall from tabled that the non-equilibrium exponents of the 
angular variables are all independent of both Tj and Tf and although the exponents do depend on d, we 
see from (I5.56|) that X^q does not. Taken literally, this would be an example of a non-universal value 
of the limit fiuctuation-dissipation ratio. 

We recall that most 'physically reasonable' systems undergoing nonequilibrium critical dynamics 
one usually finds < < 1/2, see P^l for a review. Motivated from mean-field theories of spin 
glasses, it is sometimes suggested that T^s '■= T/X^o might be interpreted as an effective temperature 
for which the fluctuation-dissipation theorem would hold. It is hard to see how in this case {Xoo may 
even become negative) such an interpretation could be maintained. 

6 The critical voter model in d dimensions 

We now study a qualitatively different type of application of local scale-invariance in the so-called voter 
model, see jHl] and references therein. The model describes the temporal evolution of conflgurations C 
of spins (Jr{t) = ±1 on a c?- dimensional hypercubic lattice Z'^. The dynamics is assumed to be given by 
a master equation 




(5.54) 




(5.55) 




(5.56) 




(6.1) 
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Here the configuration C^^^ is obtained from C by inverting tlie single spin at site r. Finally, the 
transition rates for a spin reversal a^, i— > — cr^ are given by 



1 ^ 



k=l 



(6.2) 



where the e^, 
lattice. 



d form an orthonormalized basis of unit vectors on the rf-dimensional hypercubic 



With respect to the kinetic spherical model and the XY model studied previously, the voter model is 
different since in general it does not satisfy detailed balance and therefore will not relax to an equilibrium 
state. By considering a general kinetic Ising model with a dynamics respecting the global Z2-symmetry, 
it can be shown that the transition rates ()6.2j) correspond to the critical point of the so-called linear 
voter model [7^ I88j . The non-equilibrium kinetics of the critical voter model (|6.2j) has been studied in 
detail by Dornic j2Sl- In d = 1 dimensions the model coincides with the kinetic Glauber-Ising model at 
zero temperature (which we shall revisit in appendix C) but for d > 1 these two models are different. 
In particular, it is known that the domain growth of the voter model is not driven by the minimization 
of the surface tension between the two phases ^31 1^ but which is the mechanism which drives ageing 
in simple ferromagnets lllj. 

We are interested here in the correlation functions Cr(t) = {(Jrit)ao{t)) and Crit,s) = (o-r(t)(To(s)) 
which are easily seen to satisfy the following equations of motion j2Hl I2n] 



d_ 
di 



subject to the following boundary conditions 



ArCr{t, S) 



Coit) = 1 , a(0) = 5r,0 , Cr{t,t) = Crit) 



(6.3) 



(6.4) 



where is the discrete Laplacian and where the initial magnetization '^^i'^ri^)) = 0. For the 
autocorrelation function C{t,s) = Co{t,s) of the critical voter model ()6.2|) one finds in the ageing 
regime, where t, s and t — s are all sufficiently large with y = t/s 



C{t,s) 



I arctan ^2/{y-l) 
ln(s)-i In ((y + !)/(!/ -1)) 



d-2 



-d/2+l 



(y+1 



.-d/2+1 



i{d = l 
lid = 2 
if 2 < < 4 



(6.5) 



where 7^ is the probability that a random walk in d dimensions and starting from the origin never 
returns. We are not aware of published results on i?(t, s) for 2 < d < A m. this model. Clearly, 
time-translation invariance is broken for all d > 2. 



We wish to compare ()6.5|) with the prediction eq. ()4.34|) of local scale-invariance. The case d = 1 
will be dealt with in appendix C and since for d = 2 logarithmic scaling is found, the form of local scale- 
invariance as presented here is not applicable.^ We therefore concentrate on the dimensions 2 < d < 4. 

®An extension of Schrodinger invariance to logarithmic Schrodinger invariance in analogy to logarithmic conformal 
invariance, see e.g. [44'. might be needed here. 
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From the asymptotics of C{t, s) in (|6.5|) we read off the exponents 

& = ^ - 1 , Xc = d (6.6) 
since the exponent z = 2 is known, see [Mj EE] • From eq. ()4.H4j) we expect 

fc{y) = $o,c ^'d^^ {y + i- 26)-"' = jzri[^y- - + ^y""^'] (6.7) 

in full agreement with (j6.5|) and we identify the normalization constant $o,c = ((i/27r)'^^^. Indeed, 
we see that the form (j4.33|) - which in principle is only valid asymptotically - is in fact exact in the 
critical voter model. This is not surprising in view of the underlying free-field theory. 

In conclusion, for the critical voter model with two competing steady states ageing occurs. The 
scaling form of the two-time autocorrelation function is in exact agreement with the prediction of local 
scale-invariance. This is the first time such an agreement is found for a system without detailed balance. 



7 The free random walk 

Last, but not least, we briefly consider the simplest example of a system undergoing ageing: the free 
random walk j2I] . The Langevin equation describing the time-evolution of the order parameter reads 

^ = ftW+.W (7.1) 

where a deterministic external field h{t) has been added in order to be able to compute response 
functions. The Gaussian noise t] is characterized as usual by its first two moments, see eq. Here, 
we choose the notation such that the relationship to local scale-invariance becomes evident. 

The autocorrelation and linear autoresponse functions were already calculated by Cugliandolo et al. 
They obtained, with the initial condition C(t,0) = 



R(t,s) . 



h=0 ^-^ 



6h{s) 

C{t,s) = {(j){t)(f){s)) = 2Tmin(t,s) (7.2) 

Clearly, the system undergoes ageing, since C{t,s) does not merely depend on t — s and furthermore, 
it stays forever out of equilibrium, since the fluctuation-dissipation ratio X{t, s) = 1/2 21j. 

We wish to check that the results ()7.2j) are compatible with the predictions of local scale invariance. 
For the autoresponse function, the first eq. (|7.2j) is clearly compatible with (j4. 714. 1714. 18p . with the 
exponents a = —1 and = 0. For the autocorrelation function, we expect from (j4.12p 



C{t, s)=T du 



5h{u) 



(7.3) 



where in view of the initial condition C(t, 0) = used in eq. (j7.2|) we assumed a vanishing initial 
correlator. 

In order to calculate the above derivative, we solve the Langevin equation for a given field h and 
obtain the autocorrelation function 

C{t,s-[h]) = C{t,s-[Q\)+ f [ dvdw h{v)h{w) (7.4) 

Jo Jo 
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where s; [0]) = C{t, s) is of course given by the second eq. (|7.2|) . The required second derivative of 
the autocorrelation function becomes in a field-theoretic formulation some three-point correlator (see 
eq. ()4.10p ) and it is now easy to see that 



43,(,,,,),mM;W) 



= 2e{t - u)e{s - u) (7.5) 

h=0 



5h{uf 

Inserting this into (j7.3p . the desired result for C(t, s) in eq. (j7.2p is indeed recovered. We read off the 
exponents h = —1 and Ac = 0. Of course, the exponent equalities a = b and Ac = are a necessary 
requirement for having a non- vanishing limit fluctuation-dissipation ratio Xqo = 1/2. 

In conclusion, the evidence from the two-time autocorrelation and autoresponse function of the free 
random walk is fully consistent with local scale-invariance. 



8 Conclusions and discussion 

We have analysed the ageing behaviour in systems with a non-conserved order parameter and described 
by a Langevin equation. Our main assumption was that the noiseless part of that Langevin equation 
is Galilei-invariant. Together with dynamical scaling this hypothesis fixes the dynamical exponent 
z = 2 and implies for local theories Schrodinger-invariance jlH]- There are good reasons for admitting 
such a hypothesis. For example the phase-ordering kinetics of the Glauber-Ising model in > 1 
dimensions quenched to a temperature T < Tc provides strong evidence that the scaling function of its 
space-time response function R{t, s; r, r') has the form predicted from Galilei- invariance 0^1 • However, 
since groups of local scale transformations such as the Schrodinger group are dynamical symmetries of 
noiseless differential equations only, the role of the noise in the Langevin equation or from the initial 
conditions has to be addressed. 

We have carried out such a study, for the important special case where z = 2 and the initial state is 
fully disordered. Considering the Langevin equation as the classical equation of motion of a MSR-type 
field theory, we have calculated the two-time correlation and linear response functions by studying how 
the dynamical symmetry properties of the noiseless part of that field theory are reflected in these noisy 
averages. These averages can be written in form of a perturbative expansion around the noiseless theory 
and we have shown that only a finite number of terms in these series contributes. Specifically, we have 
found: 

1. The two-time linear response function R = {(jxp) involving only quasiprimary fields is independent 
of both the thermal and the initial noise. This explains why the form ()1.7|) of scaling function 

of the linear response - previously derived from the symmetries of the noiseless theory - 
has been reproduced in many different systems with T > either exactly or with a considerable 
numerical precision [iaiIHliniEHlEnill3liailIlll3EHlinSIZZj- 

2. We obtain the reduction formula eq. (j4.9p which expresses the two-time correlation function in 
terms of certain noiseless three- and four-point response functions. For the uncorrelated initial 
conditions ()4.1H) only a single noiseless three-point response function is needed. 

3. The scaling forms of correlation and response functions are governed by the two non-trivial expo- 
nents Ac and Xr which are in general distinct from each other. Given the initial correlator ()4.11|) . 
local scale-invariance with z = 2 provides a sufficient condition for the exponent equality 

Xc = Xr (8.1) 
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4. The scaling of the two-time autocorrelator C{t, s) is described by a scahng function which in 
turn depends on a scahng function which arises in the three-point function of quasiprimary 
fields in Schrodinger-invariant theories. Depending on whether the thermal or the initial noise is 
dominant, two distinct scaling forms eqs. ()4.25|1 and ()4.32|1 are found and we have argued that 
they should describe the cases when the system is quenched to temperatures T < Tc and T = Tc, 
respectively (in the latter case only if 2; = 2 still holds after renormalization). 

5. Schrodinger-invariance by itself does not determine the form of ^E'(p). We have argued that 
the related three-point response function should be non-singular and then find the asymptotic 
behaviour for w — 00 

\ -v / - Ac ; if T < Tc , . 

$H-^ ' ' ^=1 1 + a ; if T = T. ^^-2) 

This suggests the following approximate forms, with y = t/s 

^ ^ , Mlliy + iy/iAy)]-'^^' ; if T < 

C{t,s)^{ (8.3) 

$o,cy^+"~^^/Vo^de^^c-2-2a^^_^^_20)-i-« ; if T = T, 

At least, these forms are consistent with the required asymptotic behaviour of fc{y) as ?/ — >• 00, 
see section 1. For free-field theories ()8.2j) holds for all values of w and then ()8.3|) becomes exact. 

In the past, approximate expressions for the scaling of magnetic correlation functions were derived 
from Gaussian closure procedures for kinetic 0(n)-models undergoing phase-ordering kinetics, see 
[Oj. This gives for the magnetic autocorrelation function at T = [3 El IBS! 



/c,bpt(2/) ~ 7^ 



B 



1 n 
2'~ 




lA) 



where B is Euler's beta function and 2-^1 a hypergeometric function. It is easy to see that the 
scaling function ()8.3p with T = is recovered from ()8.4|) in the n ^ 00 limit. In the notation of 
section 5.2, eq. ()8.4|) implies an exponent A'^^ = d/2. This is a typical value for a free-field theory 
for T < Tc (which indeed describes the 0(n)-model in the n — > 00 limit) but which will in general 
not hold for n < cx) and thus (|8.4p cannot be expected to represent well the behaviour ioi y = t/s 
large. ^ One might wonder whether the long-standing difficulties in arriving at scaling functions 
which cover adequately the whole range of values of y should not be related at least partially to 
having worked with dynamic variables which might turn out to be not the most basic ones of the 
model ? 

We have tested these predictions on the exact solutions of the kinetic spherical model and the 
XY model in the spin- wave approximation. In these cases, the exponent a is given by ()1.4|) . In 
order to compare the exact model results with the prediction (18. 3j) . it was necessary to carefully 
identify the quasiprimary fields of the models. For the spherical model, the natural magnetic 
order parameter and its response field could be used as quasiprimary fields. On the other hand, 
for the XY model it turned out that the coarse-grained magnetic moment is not quasiprimary but 
rather the angular variable is. 



^Indeed, to leading order in 1/n and for uncorrelated initial conditions, A'^ = d/2+{A/3)'^{d+2)^B{l+d/2, l+d/2)n~ 
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These examples underline the importance of the correct identification of the quasiprimary fields 
in a given model. 

The explicit results in the XY model also illuminate in a new way the Godreche-Luck conjecture 
on the universality of the limit fiuctuation-dissipation ratio. Further tests on the form of the 
autocorrelation function in the Glauber-Ising model in d > 2 dimensions are presently carried out 
and will be reported elsewhere [7^1 inH] • 

7. We also showed that the autocorrelator in the critical voter model in 2 < d < A dimensions and 
which does not satisfy detailed balance, again agrees with ()8.3|1 . This is the first example of a 
new domain of application of local scale-invariance. We also confirmed ()8.3|1 for the free random 
walk. 

These four confirmations, although all based on an underlying free-field theory provide further 
evidence in favour of a Galilei-invariance of the noiseless theory. 

8. The scaling of the linear response of the ID Glauber-Ising model at T = can only be explained 
through a generalization of the representations of the Lie algebra of local scale-invariance. It 
would be interesting to see whether a similar phenomenon could be found in different ID systems 
undergoing ageing at T = 0. 

Our approach has been based in an essential way on the assumption of Galilei-invariance of the noiseless 
theory. However, Mazenko jHZI recently studied phase-ordering kinetics in the context of the time- 
dependent non-linear Ginzburg-Landau equation. He carried out a second-order perturbative calculation 
around a Gaussian theory which is equivalent to the Ohta-Jasnow-Kawasaki approximation and reports 
a deviation of the two-time autoresponse function R{t, s) from the local scale-invariance prediction ()1.7|) . 
A similar difficulty had been observed before by Calabrese and Gambassi who studied non-equilibrium 
dynamics (that is T = Tc) of the 0(n) model, for both model A and model C dynamics, through MSR 
field theory ^SJEJ^j. Again, already their classical action is manifestly not Galilei- invariant. 

At face value, Mazenko's result [HZj is in disagreement with the simulational data obtained from 
the 2D and 3D Glauber-Ising model with T < T^. and based on the master equation. These do 
reproduce ()1.7p for the autoresponse function R{t, s) jl7lliniE2] as well as the extension for the spatio- 
temporal response R{t,s]r — r') Could this mean that there are subtle differences between the 
formulation of stochastic systems either through a master equation or else through a coarse-grained 
Langevin equation^^ and which affect the formal Galilei-invariance of the theory ? Alternatively, if 
under renormalization the dynamical exponent z = 2 remains constant, might the theory flow to a flxed 
point where asymptotically Galilei-invariance would hold ? 

All in all, based on a postulated extension of dynamical scaling to some local scale-invariance, we 
have reformulated the problem of flnding the scaling function of the two-time autocorrelation function 
of ageing systems as one of a discussion of the properties of certain three-point response functions of 
the noiseless theory. The evidence available at present suggests that this approach might be capable of 
shedding a new light on the issue. Finally, it might be of interest to search for extensions for dynamical 
exponents z ^ 2 and/or to study ageing systems with a conserved order parameter. Asymptotic 
information on the two-point functions in the latter case is now becoming available [10] . 



^'^At criticality, the combined effect of thermal and initial fluctuations leads for interacting theories to a non-trivial 
value of z 7^ 2 under renormalization so that our arguments are no longer directly applicable. 

^^In ID and at T = 0, the autocorrelation exponent Ac = 1 found in the Glauber-Ising model |37[ I65| differs from the 
exactly known exponent Ac — 0.6006 . . . [lOj determined in the time-dependent Landau- Ginzburg equation. 



34 



Appendix A. On Gaussian integration 

We present the details of the calculations of the magnetic two-time correlation function T and its 
associated linear response function p for the d-dimensional XY model in the spin-wave approximation. 
They are defined as 

T{t,s;r,r') = J VcpVcf) cos{(j){t,r) - ^{s,r')) exp (^-S[(p,(f)]^ 

p{t,s;r,r') = y"l)0P0 r') sin r) - r')) exp 0]) (Al) 

where S[(f),(j)] is the free- field Martin-Siggia-Rose action which is Gaussian. 

A.l The correlation function r{t, s;r,r') 
Using de Moivre's identities 

T{t,s-r,r') = ^[(expi(0(t,r)-0(s,rO)) + (expi(0(s,r')-0(t,r)))] 



exp i J du dR J{u, R)(j){u, R)^ + ^exp — i J du dR J{u, R)(j){u, R) 



(A2) 



where J{u, R) is given by 

J{u, R) = 6{u - t)6{R -R)- 6{u - s)6{R - r') (A3) 
A generally valid result from free-field theory reads, see 



expi JdudRJ{u,R)(p{u,R)^ = exp j dudu' dRdR' J{u, R) {^{u, R)(l){u' , R')) J{u , R') 

(A4) 

which is one of the various forms of writing Wick's theorem. Using the explicit form of the current in 
(IA3jl . equation ()5.33|1 follows immediately. 

A. 2 The response function p{t, s;r,r') 

We now focus on the magnetic response function p{t, s; r, r'). We decompose the MSR action 

5[0,0,O] = 5o[0,0] + s[0,0] (A5) 

where 

5o[0,0] = Jdudr^ ^- 
s[<^, (/)] = — j dudu' dr dr' (j){u,r)K{u,u']r — r')(j){u' ,r') (A6) 



and 



k{u, u'] r-r)= T6{u - u')6{r - r') + ^6{u)6{u')a{r - r') (A7) 
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(r-r'Y 



A(t - s) 



The noiseless two-point functions are 

Ro{t,s-r,r') = r)0(s, r'))o = e(t - s)[47r(t - s)]-'^/^ ^xp 

Co{t,s-r,r') = (0(t,r)0(s,r'))o = O 

Co(t,s;r,r') = r)0(s, r'))o = (A8) 

Below, we shall need the equal-time response function s; r, t')|j^^. To give to this quantity a value, 
one may discretize the Langevin equation. This may be done according to several different schemes, 
see jnniinS- Here, we shall use the Ito prescription which amounts to 

Ro{t,t-r,r')=Q (A9) 

The magnetic response function reads from ()A1|) 

p(t, s- r, r') = r') sin(0(t, r) - r')) exp (AlO) 



where ()o is the average with the exp ^— S'o[0, 0]^ . Expanding the sine and using the Newton's binomial 
identity, we have 

oo 2fc+l 

p(t,s;r,r') = ^ ^ pfc,p(t, s; r, r') 

fc=0 p=0 

p,,,{t,s;r,r') = ^,(2^ + 1 - p)\ V^'^ r)0^^-+^-^(., r') exp -.[0, 0])^ (All) 

so that Pk,p{t, s] r, r') is given by a 2A; -|- 2-point function, containing only one (p contribution and 2k + 1 
fields 0. The Bargman superselection rule ()A8|) implies that the only contractions that will lead to 
non- vanishing averages come from the k-th order in the expansion of the exponential, so that we have 

p.,,(t,.;r,rO = ^,^J~^]'^'^^,^, J f[djdf (0(., r')0^(a;, t)02^+^-^(., r')0(j>U (A12) 

where for clarity the notation j, (respectively j') stands for {uj,rj)^ (respectively {upr'j)) and where 
integrals run over this set of 2k variables. 

Wick's theorem states that the integrand decomposes into sums of products of two-point functions. 
In order for a contraction not to vanish, the field 0(s, r') must contract with one of the p fields 0(t, r), 
which leads to 

Pk,Q{t,s;r,r') = 

x(^r-\t,r)<p''^'-^is,r'mj)Kij,j')^{j')^ (A13) 
Summing over p, the response function p(t, s; r, r') reads 

p(t,3;r,r')=i?o(t,3;r,r')5^-TT^ / {m,r) - <j>{s,r')f'^{j)K{j, f)^{f)) (A14) 

A:=0 ^ ■ j=l " 
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At this stage, it will be interesting to give another equivalent expression of the magnetic correlation 
function r(t, s; r, r'). Using the same strategy as before, but now expanding the cosine we find 

oo 

r(t,s;r,r') = ^ 7fc(t, s; r, r') 

fc=0 

lk{t,s-r,r') = tl)!(^(0(t,^)_0(s,r'))2'=exp-s[0,0])^ (A15) 

For the same reason as before, because of mass conservation, only the expansion of order k of the 
exponential term will contribute and lead to non-vanishing contractions, such that 

k 

7k{t,s;r,r') = |ndjd/ r) - r'))^^0(j>(j,/)0(/))^ (A16) 

Comparing ()A14|) with ()A16|) . we thus have 

p{t, s; r, r') = Ro{t, s; r, r')r{t, s; r, r') (A17) 
This is the main result of this appendix and gives eq. ()5.38|) in the text. 



Appendix B. Scaling form of a special four-point function 



We study the scaling of the four-point function 

:= R^o^ {t, s, 0, 0; 0, 0, R, R') = ( O)0(s, 0)0(0, il)0(O, R 



(Bl) 

of a theory in MSR formulation of which the noiseless part is Schrodinger- invariant. The field (j) is 
assumed to be quasiprimary with scaling dimension x and mass Ai and the response field (j) should also 
be quasiprimary with scaling dimension x and mass Ai = —M.. The covariance conditions on TZ are 
the following (we use v{t) = 0) 

(^tdt + sd, + ^RdR+^Rdji> + {x + x)^J^ = 

(^t^dt + s^ds + {t + s)x-^(^R^ + R'^^^J^ = 



(B2) 



We shall not discuss here the general solution of these equations. For our purposes it is enough to 
observe that if we decompose in the following symmetrized way 



= g{t, R)g{s, R) + g{t, R)g{s, r) 

which for a free-field theory would follow from Wick's theorem, then 



e?(t,i?) = 6;ot~"""exp 
produces a solution to both covariance conditions. 



MR^ 



(B3) 



(B4) 
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We apply this result to the preparation part Cpr of the autocorrelation function, with the intention 
to use the result in the spin-wave approximation of the XY model. From eq. ()4.9|) we have 



pr 



- / dRdR' a{R - R')R^Q\t, s, 0, 0; 0, 0, R, R) 



(2vr; 

oo 



(47r 



d 

dz 

]-d/2 



dRdRdq ^ ' „ exp 



; ,„ M [R^ R 



/2 



dq 



(B5) 
(B6) 
(B7) 
(B8) 



d/2 - 1 

where in going to ()B6|) we used the initial condition (|5.27|) and the result ()B3IB4|) from above, next in 
going to ()B7|) we speciahzed to a free Gaussian field (where x = x = d) oi mass Ai = 1/2 and the last 
step we also assumed d > 2. 



Eq. ()B8|) provides the preparation term, for any initial temperature T,, as required for the analysis 
of the autocorrelation function in the XY model in section 5.2 



Appendix C. On local scale- invar iance in the ID Glauber-Ising 
model 

C.l Two-point functions in the ID Glauber-Ising model 

In the text we have seen that local scale-invariance implies the following form of the two-time autore- 
sponse function 



R(t,s)=ro{t-s)-'-''\^-j (CI) 

In spite of a nice agreement with a large variety of models, this expression is not verified for the ID 
Ising model with Glauber dynamics at T = 0. 

The ID Ising model is described by spins cij = ±1 and the Hamiltonian 

n = -Y^cx,a,+, (C2) 

i=l 

The exactly solvable Glauber dynamics [SEj niay be given through the heat-bath rule, which gives the 
probability of finding the spin variables crj(t + 1) in terms of those at time t 

P {ai{t + 1) = ±1) = i [1 ± tanh (/3 (a,_i(t) + (Ji+i{t) + h,{t)))] (C3) 

where = 1/T is the inverse temperature and hi{t) the external magnetic field. Ageing occurs in this 
model at T = 0. In the long-time scaling limit, the two-time autocorrelation and autoresponse functions 
are jSa EHl EHl 

1 



Shi(s) 



h=o 7T^/2s{t-s 



(C4) 



C{t,s) = {ai{t)ai{s)) = ^ arctan ^ ^ (C5) 
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While these results were obtained first for a fully disordered initial state, they remain true for long- 
ranged initial conditions ((Tr(0)(To(0)) ~ r'" with > ||51i- The case i/ = corresponds to the case 
of an initial magnetization rrtQ. Then the connected part of C(t, s) as well as R{t, s) are multiplied by 
1 — ml In any case, the forms of the scaling functions fcniv) are unchanged. 

Although these two-point functions clearly display dynamical scaling, it is evident that the scaling 
form of R{t,s) from ()U4|) is incompatible with the form suggested in ()C1|) . Local scale- invariance as 
developed in the text does not hold in the ID Glauber-Ising model. 



C.2 Generalized realization of the ageing algebra 



We now show how Schrodinger invariance can be generalized such that the exact response function ()C4|) 
can be reproduced. Obviously, time-translation invariance is broken in ageing systems. Therefore, as 
already pointed out in 1^302]; the dynamical symmetry cannot be the Schrodinger Lie algebra sci)i 
which contains the time-translation generator X_i = —dt, but a subalgebra without this generator 
might be acceptable. We consider the algebra 



{Xo,Xi,F_i/2,Fi/2,Mo} 



(C6) 



and keeping the commutation relations (j2.7|) we now look for a more general realization of agc^. In this 
way, we write the generators as {Sq^i, T±i/2, Mq}. These must be of the form 



-tdt rdr — — 

2 2 



trdr — xt — g{t) — -^'^^ 



Mo 



-dr 

-tdr - Mr 
-M 



(C7) 



where g = g{t) is to be found. The only commutator of aQti constraining g is 



-1) --oj 







(C8) 
(C9) 



which leads to 

tdtg - g 

with the solution g{t) = Kt, with K some constant. From ()C7|) . the dynamical exponent z = 2. If we 
were to require in addition = 2So (and thereby go from oge^ back to scf)i) we would recover 

K = 0. Now, a quasiprimary field of age^ will be characterized by a triplett {x, K, Ai). 

We can now generalize local scale-invariance by requiring that the autoresponse function R{t, s) 
formed from a quasiprimary field (p and its associated quasiprimary response field (p to transform 
covariantly under the generators Sq and Si. It is a solution of the system of linear partial differential 
equations 



t^t + - + sds + - 



R{t,s) 



t'^dt + {K + x)t + s'^ds + (K + x)s R{t,s) 








(CIO) 
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where x and K refer to the response field 0. Solving the system ()C10|) gives as final result 

R{t,s) = s-'-'^fnit/s) 

fdy) = (y - 1)-''^ (Cll) 

where the three independent non-equilibrium exponents a, A, A/j are 

X + x 



A = a+K+K 

\b = 2x + 2K (C12) 

In contrast with the previous realization of agc^, K,K is possible and then a and A differ from 
each other. 

Comparison with the exact result ()C4|) of the ID Glauber-Ising model now gives complete agreement 
and we identify the exponents 

a = , A = -^ , Xn = l (C13) 

Of course, the values of a and Xr have been obtained before but A seems to be a new exponent. 

At present, it must remain open whether the unusual properties of the ID Glauber-Ising model are 
related to the fact that = and therefore the critical and low-temperature properties might have 
become mixed. 

Also, it remains to be seen whether the form of the autocorrelation function can be understood from 
the generalized realization of the ageing algebra agc^. We hope to come back to this elsewhere. 
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